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1.0  INTRODUCTION 


The  solar  wind  (a  charge  neutral  plasma  comprised  mostly  of  protons  and 
electrons)  is  the  ultimate  energy  source  for  many  phenomena  In  the  earth's 
upper  atmosphere  and  ionosphere--especially  at  high  latitudes.  The  solar 
wind  also  is  entirely  responsible  for  the  formation  and  maintenance  of  the 
earth's  magnetosphere.  These  regions  of  the  near  earth  environment  are 
strongly  coupled  and  in  turn  all  act  to  influence  space  and  military  hard- 
ware systems  in  many  ways.  Examples  include:  the  influence  of  the  ionosphere 

on  radio  communications  and  the  uses  of  over-the-horizon  radar;  the  influence 
of  the  neutral  atmosphere  on  satellite  trajectories  and  lifetimes;  the  effects 
of  trapped  and  cosmic  radiation  on  man  and  components  in  space;  and  the 
influence  of  space  currents  on  surface  systems  vulnerable  to  magnetic  field 
fluctuations.  In  order  to  optimize  the  use  of  these  hardware  systems,  it  is 
necessary  to  understand  quantitatively  these  environmental  influences.  For 
the  real  time  optimization  of  several  systems  it  is  necessary  to  predict  the 
behavior  of  the  environment  in  which  these  systems  operate.  At  present, 
most  Air  Force  efforts  In  this  area  are  limited  to  the  real  time  "specification" 
of  environmental  features.  Air  Force/Global  Weather  Central  (AFGWC)  routinely 
receives  real  time  data  from  several  sources  which  are  used  for  specification 
purposes.  However,  to  date  It  has  not  been  possible  to  offer  a bona  fide 
predictive  capability  (except  for  the  forecast  of  a few  space  averaged  bulk 
parameters). 

There  are  several  reasons  for  the  present  lack  of  predictive  capability.  The 
several  regions  of  the  near  earth  environment  and  their  coupling  are  not  cur- 
rently adequately  understood.  The  energy  sources  that  produce  structures  and 
variability  in  the  near  earth  environment  must  be  monitored.  At  present,  this 
Is  not  done  routinely.  Quantitative  models  of  the  environmental  parameters  In 
question  and  the  relations  to  their  energy  sources  must  be  developed  with 
associated  fast  computer  codes. 

In  this  contract,  the  response  of  the  earth's  magnetosphere  and  upper 
atmosphere  to  the  solar  wind  has  been  studied.  Several  goals  were  stated  and 
met.  The  overall  purpose  of  the  work  was  to  obtain  enough  understanding  of 
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some  environmental  features  to  develop  quantitative  models  of  their  behavior, 
w This  was  done  with  the  intent  of  examining  the  possibility  of  predicting  solar 
wind  influences  on  the  near  earth  environment. 

Four  general  problem  areas  were  isolated  with  the  above  goal  in  mind. 

1.  Work  prior  to  this  contract  was  done  on  energy  sources  for  the  upner 
atmosphere  (Moe,  1971;  Olson,  1971).  It  was  found  in  particular  that 
at  high  latitudes  charged  particles  from  the  solar  wind  precipitate 
directly  into  the  ionosphere  and  upper  atmosphere  near  noon  (Frank,  1971, 
and  Heikkila,  1971).  This  precipitation  is  controlled  by  the  magneto- 
spheric  magnetic  field.  The  "field  lines"  along  which  these  particles 
move  form  the  region  of  the  magnetosphere  known  as  the  dayside  cusp  or 
cleft.  The  physics  of  the  dayside  cusp  was  studied  during  the  first 
year  of  the  contract  as  were  the  correlations  discussed  above.  Data  on 
precipitating  particles  were  used  to  describe  the  extent  and  location  of 
the  cusps.  The  amount  of  corpuscular  energy  deposited  through  the  cusps 
was  also  determined  as  it  is  an  energy  source  for  both  the  upper  atmos- 
phere and  the  Ionosphere  and  thus  an  important  input  parameter  for 
quantitative  models  of  those  regions.  The  currents  in  the  cusp  region 
and  their  associated  fields  were  also  estimated.  The  cusp  work  is 
discussed  in  Section  2.0. 

2.  Any  attempt  to  predict  atmospheric  and/or  ionospheric  behavior  must 
consider  the  coupling  of  both  regions  to  the  magnetosphere.  This  is 
necessary  because  solar  wind  variability  must  be  monitored  in  advance 
of  the  time  for  which  near  earth  environmental  conditions  are  being 
predicted.  The  solar  wind  interacts  with  the  earth's  magnetic  field 
to  form  the  magnetosphere.  The  magnetos pheric  magnetic  field  in 
particular  exerts  considerable  Influence  on  the  near  earth  environment, 
especially  at  high  latitudes.  Therefore,  an  attempt  was  made  to  more 
accurately  model  the  total  magnetos pheric  magnetic  field.  By.  By  has 
as  its  main  sources  the  earth's  main  (Internal)  field,  the  magnetopause 
field  (from  currents  formed  by  the  deflection  of  the  shocked  solar  wind 
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electrons  and  protons  In  opposite  directions  by  By),  the  tail  field 
(from  particles  that  drift  across  the  tail),  and  the  ring  current 
magnetic  field  (produced  by  particles  trapped  in  By  and  drifting  around 
the  earth).  Earlier  quantitative  models  of  By  had  neglected  the  ring 
and  tail  currents  since  it  is  very  difficult  to  modlfK  magnetic  field 
in  the  region  of  its  source  currents.  (Scalar  potentials  cannot  be 
used  in  such  regions.)  A method  was  developed  to  overcome  these  problems. 
The  result  was  the  first  accurate  quantitative  model  of  By  that  included 
all  four  of  its  major  sources.  This  model  is  discussed  in  Section  3.0. 
This  work  was  performed  during  the  first  and  second  years  of  the  contract 
period. 

3.  The  studies  on  the  dayside  cusps  (or  clefts)  confirmed  what  many  earlier 
investigations  had  suggested--that  at  high  latitudes  the  upper  atmosphere 
is  very  structured  and  charged  particles  are  an  important  energy  source 
in  that  region.  These  conclusions  led  to  an  attempt  to  construct  a 
quantitative  global  model  of  the  upper  atmosphere  that  would  for  the 
first  time  include  these  high  latitude  corpuscular  energy  sources.  This 
model  and  the  model  of  By  are  semiempi rical . That  is,  they  use  physical 
understanding  whenever  it  exists,  but  rely  on  observational  data  in  the 
many  areas  where  the  physics  is  currently  not  well  known.  Work  on 
modeling  the  density  of  the  upper  atmosphere  was  performed  during  the 
second  and  third  years  of  the  contract.  The  atmospheric  density  model  is 
discussed  in  Section  4.0. 

4.  An  investigation  of  the  several  reported  correlations  between  solar 
wind  apd  near  earth  environmental  parameters  was  carried  out.  Its 
purpose  was  twofold.  First,  an  attempt  was  made  to  isolate  highly 
correlated  data.  This  would  ultimately  suggest  a means  by  which  some 
near  earth  environmental  parameters  could  be  predicted  directly  from 
monitored  solar  wind  data.  Such  correlations  also  suggest  that  the 
direct  monitoring  of  solar  wind  data  might  yield  Indices  that  are  much 
more  accurate  and  useful  than  the  temporally  and  spatially  averaged 
indices  such  as  Kp,  A£,  and  DSy.  Second,  such  correlations  must 
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ultimately  be  used  to  guide  the  development  of  quantitative  environmental 
models.  Such  models  are  only  useful  to  the  extent  that  their  input 
parameters  are  available  (e.g.,  from  a monitoring  satellite)  and  their 
output  parameters  accurately  describe  environmental  features  (e.g., 
atmosphere  density).  These  correlations  are  discussed  in  Section  5.0. 

In  addition  to  the  four  general  problems  discussed  above,  there  were  several 
other  areas  of  direct  or  indirect  results  associated  with  performance  on  this 
contract.  Both  models  were  developed  with  the  user's  needs  in  mind.  This  is 
discussed  in  the  section  on  modeling  philosophy  in  Section  6.0.  Other  work 
that  has  grown  out  of  this  contract  work  is  also  discussed  in  Section  6.0. 

A summary  of  this  work  and  conclusions  regarding  it  are  also  presented  in 
Section  6.0.  Appendicies  to  this  report  appear  in  Section  7.0.  They  include 
listings  of  the  quantitative  magnetic  field  and  neutral  atmospheric  density 
models.  Copies  of  Fortran  decks  for  these  models  are  available  upon  written 
request  from  McDonnell  Douglas  Astronautics  Company.  (Contact  Dr.  W.  P.  Olson 
or  Ur.  K.  A.  Pfitzer  for  details:  (714)  896-4368.)  The  technical  references 
are  given  in  Section  8.0  while  the  list  of  publications,  reDorts,  and  presen- 
tations resulting  from  this  contract  is  given  in  Section  9.0. 
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2.0  DESCRIPTION  OF  THE  PROPERTIES  OF  THE  DAYSIDE  CUSPS 
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l 
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' 2.1  Observations  of  Cusp  Parameters 

Work  on  the  physics  of  the  dayside  cusps  was  Initiated  with  the  following 
information  at  hand.  The  cusp  at  its  Intersection  with  the  Ionosphere  extends 
a few  degrees  In  latitude  and  8-12  hours  In  local  time  with  its  center  on  the 
noon  magnetic  meridian.  In  the  outer  regions  of  the  magnetosphere  the  cusp 
has  flared  so  that  tt  extends  several  earth  radii  near  local  noon.  The 
position  of  the  cusp  is  determined  by  the  axis  of  the  earth's  dipole  maanetic 
field  and  the  direction  of  the  solar  wind.  Since  the  magnetic  and  rotation 
axes  of  the  earth  are  inclined  by  about  11.7°,  the  location  of  the  dayside 
cusp  in  geographic  coordinates  Is  continuously  changing. 

Particle  precipitation  through  the  cusps  yield  a height  Integrated  energy  flux 
of  -v.8  ergs  (cm  ster  sec)'  during  quiet  magnetic  periods,  Olson  (1971). 

(This  number  was  obtained  from  the  analysis  of  the  published  pitch  angle  and 
particle  spectra  of  Frank  (1971)  and  Helkkila  (1972).)  During  disturbed 
periods  the  energy  flux  can  be  much  larger.  The  protons  carry  most  of  the 
energy.  Their  peak  energy  transfer  to  the  atmosphere  occurs  at  about  170  km 
while  the  peak  electron  energy  transfer  to  the  atmosphere  is  at  about  200  km 
(Olson,  1971).  The  shape  of  the  Intersection  of  the  dayside  cusp  with  the 
atmosphere  Is  observed  to  be  (roughly)  a semi-circular  region  centered  on 
the  magnetic  pole  and  extending  about  180°  in  longitude  with  the  central 
longitude  at  magnetic  noon  (see  Figure  2-1).  In  magnetic  coordinates  the  shape 
of  the  cusp  Is  almost  constant  during  quiet  magnetic  conditions.  (It  exhibits 
some  dependence  on  the  angle  between  the  solar  wind  direction  and  the  dipole 
axis.)  In  geographic  coordinates,  however,  the  cusp  position  is  constantly 
moving.  Thus,  the  amount  of  time  a given  point  will  be  within  the  cusp 
region  and  the  extent  to  which  It  will  be  Influenced  by  charged  particle  fluxes 
depends  on  Its  geographic  location,  both  its  latitude  and  longitude.  The 
center  of  the  cusp  Is  at  about  75°  magnetic  lal!tude  during  quiet  times.  The 
seasonal  variation  in  the  cusp  location  Is  quite  small  (In  magnetic  coordinates). 
The  region  of  maximum  exposure  of  the  atmosphere  and  Ionosphere  to  the  cusp 
particles  is  at  the  dipole  longitude  (-69°  West  In  geographic  coordinates)  and 
about  15°  below  the  dipole  axis  for  the  northern  hemisphere. 
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FIGURE  2-1 


The  geographical  region  that  spends  the  most  time  under  the  cusp  (and  therefore 
receives  the  most  energy)  is  near  Hudson  Bay.  This  is  the  location  of  the 
"fixed"  atmospheric  density  bulge  reported  by  Jacobs  (1967).  His  work  was 
based  on  satellite  drag  studies.  Other  better  and  more  recent  atmosnheric 
density  data  also  suggest  heating  and  Increased  density  In  this  region. 

Since  the  energy  input  to  the  atmosphere  via  charged  particles  through  the 
dayside  cusps  is  under  magnetic  control  it  is  possible  to  distinguish  this 
energy  source  from  the  solar  UV  heating  that  occurs  on  the  dayside  of  the 
earth.  Near  winter  solstice  the  cusp  deposits  Its  energy  on  the  nightside 
of  the  earth.  In  Figure  2-2  the  zenith  angle  of  the  cusp  is  shown  at  different 
times  of  the  year.  A zenith  angle  of  greater  than  90°  indicates  that  the 
center  of  the  cusp  (at  Its  Intersection  with  the  atmosphere)  Is  on  the  nightside 
of  the  earth.  This  geometric  feature  makes  the  direct  testing  of  the  extent 
and  position  of  the  dayside  cusp  energy  source  possible.  This  has  been  done 
by  comparing  the  cusp  location  (given  by  this  model)  with  airglow  data  obtained 
over  the  northern  polar  cap  in  winter  (see  Section  5.0). 

The  model  discussed  above  permits  comparison  of  data  with  cusp  location  and 
extent.  In  order  to  answer  more  detailed  questions  about  the  formation  or 
maintenance  of  the  cusp,  a physical  examination  of  the  particles  that  move 
through  the  cusp  was  necessary. 

2.2  The  Physics  of  the  Dayside  Cusp 

The  presence  of  the  cusped  geometry  In  place  of  the  neutral  point  geometry 
obtained  from  the  older  models  had  to  be  accounted  for.  For  such  details 
as  the  cusp  geometry,  the  behavior  of  Individual  particles  had  to  be  considered. 
In  particular,  we  found  that  whenever  strong  gradients  In  BT  exist  just  Inside 
the  magnetopause,  particles  can  drift  Into  the  magnetosphere.  B exhibits  this 
geometry  In  the  cusp  region  and  In  the  vicinity  of  the  equator  along  the  tall 
of  the  magnetopause.  Particle  entry  In  the  tall  of  the  magnetopause  by  this 
mechanism  Is  believed  responsible  for  the  cross  tall  currents  (Olson,  1972) 
and  much  of  the  structure  of  the  plasma  sheet  (Pfltzer  and  Olson,  1973).  Thus, 
even  If  a model  magnetosphere  exhibits  a closed  field  line  geometry  (all  of 
the  field  lines  confined  within  the  magnetosphere).  It  Is  certain  that  the 
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DAILY  AND  SEASONAL  VARIATION  IN  ZENITH  ANGLE  OF 
CENTER  OF  DAYSIDE  CUSP  AT  ATMOSPHERIC  HEIGHTS 


FIGURE  2-2 


drift  of  individual  particles  will  make  the  magnetosphere  "open"  to  charged 
particles..  That  is,  some  of  the  shocked  solar  wind  plasma  (magnetosheath 
plasma)  will  penetrate  the  magnetosphere.  This  is  especially  true  in  the 
dayside  cusp  regions  and  the  equator  of  the  t’11.  We  concluded  that  the 
magnetosphere  cannot  be  adequately  described  ir  terms  of  pressure  balance  and 
these  individual  particle  motions  must  be  considered. 

The  magnetosheath  particles  Initially  penetrate  the  magnetosohere  In  a large 
region  (around  the  neutral  points  In  the  old  models).  Once  they  enter  the 
magnetosphere  they  obey  the  adiabatic  Invariants  and  move  along  field  lines 
with  most  of  them  precipitating  Into  the  atmosphere  or  being  mirrored  back 
into  the  magnetosheath.  Within  the  cusp  the  motions  of  the  particles  are 
diamagnetic,  the  particles  circle  the  field  lines  such  that  the  currents 

4 

they  produce  tends  to  weaken  the  total  field,  By.  The  currents  of  the 
individual  particles  tend  to  cancel  each  other  except  at  the  cusp  boundary. 

-4 

The  cusp  currents,  J , then  flow  In  loops  only  on  the  cusp  boundary.  (This 

t 

cusp  current  system  is  diamagnetic,  the  particles  that  make  up  the  current 
are  "stationary",  they  move  only  around  and  up  and  down  field  lines  - not 

•4  -4 

around  the  cusp.)  J produces  a magnetic  field  that  weakens  B within  the 

C 

cusp.  In  terms  of  field  components  measured  at  the  earth's  surface  (Z  or  V), 

4 

the  field  from  Is  in  the  positive  vertical  direction.  This  suggests  that 

C 4- 

the  polar  surface  magnetic  field  from  J should  exhibit  a dally  variation  In  Z 
or  V as  the  earth  turns  under  the  dayside  cusp  currents. 

Near  the  magnetopause  Jc  Is  quite  weak  because  many  of  the  particles  move 
close  to  the  direction  of  the  field  lines.  As  the  particles  move  down  the 
field  lines  (toward  the  atmosphere)  they  increase  their  pitch  angles  (the 

4 -4 

angle  between  the  particle  velocity  vector  V,  and  B)  until  they  "mirror" 
although  many  of  them  precipitate  Into  the  atmosphere.  The  increased  average 

-4 

pitch  angle  at  lower  altitudes  Increases  the  magnitude  of  J.  The  increased 

4 4 

particle  density  also  Increases  J.  However,  J does  not  Increase  so  rapidly 
as  to  cancel  B within  the  cusp  (like  the  magnetopause  currents  do  when  pressure 
balance  formalism  1$  used).  Thus,  although  B Is  somewhat  weakened  (especially 
in  the  mid-altitude  cusp)  B Is  finite  throughout  the  cusp. 
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The  determination  of  Jc  was  guided  by  an  empirical  knowledge  of  n,  V , V , A, 
and  B which  respectively  are  the  particle  number  density,  particle  velocity 
components  parallel  and  perpendicular  to  B,  the  cross  sectional  area  of  the 
cusp  and  the  magnitude  of  the  magnetic  field. 

In  a region  of  the  cusp  when  no  mirroring  occurs,  conservation  of  total  charge 
yields 


n V A = Constant  (2-1) 

while  conservation  of  magnetic  flux  gives 

AB  = Constant.  (2-2) 

It  is  assumed  that  A is  perpendicular  to  B and  that  the  edges  of  the  cusp  are 
parallel  to  B.  Near  the  magnetopause  the  following  values  are  observed  for 
several  of  the  variables  listed  above.  The  subscript  m indicates  a variable 
near  the  magnetopause  end  of  the  dayside  cusp  while  the  subscript  g indicates 
the  value  of  a variable  at  the  atmospheric  intersection  with  the  cusp.  From 
observations 


n 2 1 cm"3,  Vm  = 400  km  sec"1  = V 
m m 

(V  is  the  total  velocity  - since  all  processes  under  consideration  are 
conservative,  V is  constant) 

A = 4.1018cm2,  B_  ■ 10y,  while  near  the  earth 
m m 

Bg  ^ 104,  and  A^  i .05R2  (where  R^  is  the  radius  of  the  earth). 

With  the  known  values  of  several  variables  the  constants  in  equations  (1)  and 
(2)  can  be  evaluated  near  the  magnetopause. 

nm  Vm  \ “ K6*1q26  sec  l (2_3) 

8 A 2?  4.1*1014  gauss  cm2 
mm 


Also 


(2-4) 


* ' Finally  the  first  Invariant  constant  is 

k = V2/B  = 10®  cm  sec"^  gauss'^.  (2-5) 

With  (1)  and  (2)  and  the  known  values  of  several  variables  it  is  possible 
to  determine  the  strength  of  Jc  as  a function  of  altitude  in  the  cusp.  The 
number  of  particles  that  pass  through  a unit  length,  £,  along  the  cusd 
surface  is  given  by: 


n^  = £ R • 2irr*n 

where  R is  the  particle  gyroradius  and  r is  the  radius  of  the  cusp.  The 

p 

average  velocity  of  the  particles  is  V Thus,  the  total  flux  of  charges 
per  second  past  £ Is 

x2  1 m Vx 

0 = N£  . q ^ Since  R = 


where  m Is  the  particle  mass  - and  the  (constant)  first  adiabatic  Invariant 

2 

is  defined  by  p = m Vj_/2B,  the  current  per  unit  length  Jc  * j/£  becomes 


i ==  *20. 

Jc  ir  * 


(2-6) 


Since  p Is  constant  j depends  only  on  n.  The  change  In  n with  altitude  in 
the  cusp  can  then  be  determined  (1)  provided  the  dependence  of  V„  on  altitude 
is  known.  Since  V£/B  Is  constant  (a  form  of  the  first  adiabatic  Invariant) 
and  V Is  constant,  we  find  V2  = V2  - V2  = V2  - kB  where  k = V2/B,  thus 
Vn  = (V2-kB)^2  and  from  (1)  n(V2-kB)^2  A = 1.6  x 1026.  A varies  as  the 
inverse  of  B according  to  equation  (2).  However,  It  Is  realized  that  equation 
(2)  provides  only  a guideline  in  a region  where  currents  are  flowing  and  that 
the  product  AB  Is  not  precisely  constant.  Combining  these  results  yields 


n 


3.8  x 1011 


or,  when  the  constants  are  used 


n » 3.1018,5B/(1.6-3B)1/2 


(2-7) 
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when  B • .3  n = 1.8  x 1017’5 
when  B = .0001  n ^ 2 x 1014’5 

Thus  n varies  by  almost  3 orders  of  magnitude  in  this  simple  analysis. 

Actually,  the  currents  on  the  cusp  boundary  and  the  distribution  of  pitch 

angles  at  any  given  altitude  within  the  cusp  tend  to  uiminish  this  range  of 

densities.  Although  the  solar  wind  has  a directed  velocity  component,  the 

particles  do  not  all  enter  the  cusp  moving  along  field  lines  as  has  been 

assumed  here.  Further,  if  the  particle  distribution  was  isotropic,  there 

would  be  no  Increase  in  n toward  the  base  of  the  cusp.  An  intermediate  case 

where  the  density  increases  by  a factor  of  10  or  15  from  the  ton  to  the  base 

of  the  cusp  is  not  unreasonable.  Similarly  the  strength  of  J should  then 

-*•  c 

increase  by  the  same  amouni . The  value  of  B near  the  currents  (at  a point 
*► 

close  enough  to  J that  the  currents  appear  to  be  planar)  is  linearly  pro- 
portional  to  J.  Thus,  if  B is  7y  near  the  magnetopause  it  will  be  about  100y 
in  the  ionosphere.  At  points  more  distant  than  this  B will  be  somewhat 
smaller.  At  the  base  of  the  cusp  then  B will  have  a component  of  about  lOOy 
opposite  to  the  dipole  contribution.  Thus,  a ground  based  magnetometer  should 
measure  a dally  variation  on  the  order  of  100y  at  points  that  pass  under  the 
dayside  cusp.  Such  behavior  in  the  polar  surface  magnetic  field  has  long  been 
observed.  The  largest  dally  variations  occur  at  the  cusp  latitude  - about  75° 
magnetic  latitude.  Also,  the  daily  variation  has  the  proper  phase.  The 
vertical  component  (positive  out  from  the  earth)  is  largest  near  noon. 
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3.0  QUANTITATIVE  MODEL  OF  THE  TOTAL  MAGNETOSPHERIC  MAGNETIC  FIELD 


In  order  to  use  solar  wind  data  to  predict  atmospheric  and  ionospheric 
behavior  it  is  necessary  to  have  a quantitative  model  of  the  total  magneto- 

•f 

spheric  magnetic  field,  B,.  Charged  particles  precipitating  through  the 
dayside  cusps  and  solar  cosmic  rays  are  controlled  by  By.  The  model  of  By, 
therefore,  should  be  capable  of  accepting  as  input  descriptors,  several  solar 
wind  parameters  capable  of  being  monitored  by  satellite  borne  sensors.  Out- 
put from  the  model  can  then  be  used  to  predict  the  location  and  extent  of  the 
dayside  cusp  regions  and  the  variation  in  those  other  field  locations  that 
control  particle  precipitation  from  the  tail  and  middle  maqnetosphere. 

Existing  magnetic  field  models  have  proven  inaccurate  for  many  applications. 

A study  of  the  observations  of  the  magnetos pheric  magnetic  field  was  there- 
fore made  in  order  to  determine  those  regions  of  the  magnetosphere  where 
existing  models  were  grossly  In  error.  Two  important  conclusions  were 
reached.  First,  on  average,  in  the  tail  of  the  magnetosphere,  field  lines 
from  the  two  lobe  regions  are  connected  to  the  equatorial  region  of  the  tail. 
In  earlier  models  an  infinitesimally  thin  current  sheet  flowed  across  the 
equatorial  region  of  the  tail  and  prohibited  the  return  of  field  lines 
through  the  equator.  Second,  observations  of  the  magnetic  field  in  the 
inner  magnetosphere  (In  the  range  of  geocentric  distances  from  2 to  6 earth 
radii)  indicate  that  even  during  quite  magnetic  times  the  magnetic  field  is 
depressed.  This  Indicates  that  a distributed  ring-like  current  flows  through 
the  inner  magnetosphere  at  all  times.  The  observed  field  then  Is  better 
represented  by  a simple  dipole  main  field  with  no  contributions  from  the 
magnetosphere  than  by  some  of  the  older  models  which  incorrectly  represented 
this  region  as  being  one  of  augmented  magnetic  field.  These  two  features  are 
of  major  importance  for  many  magnetospheric  phenomena.  The  weak  field  region 
in  the  Inner  magnetosphere  will  cause  field  lines  at  a given  equatorial 
location  to  Intersect  the  earth  at  lower  latitudes  than  given  by  previous 
models.  Because  of  this,  the  poleward  extent  of  the  trapping  region  for 
adiabatic  particles  and  the  cutoff  for  solar  cosmic  ray  particles  should  be 
at  lower  latitudes  than  those  given  by  previous  models.  The  return  of  field 
lines  through  equatorial  regions  of  the  tail  dictates  that  charged  particles 
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cannot  gain  access  to  the  inner  magnetosphere  via  that  region  of  the  tail  as 
has  been  assumed  many  times  in  work  on  solar  cosmic  rays. 

In  order  to  use  this  information  on  the  observed  magnetospheric  magnetic  field 
to  develop  a magnetic  field  model  it  was  first  necessary  to  try  to  understand 
the  currents  responsible  for  the  magnetic  fields.  Work  on  this  aspect  of  the 
problem  resulted  in  a model  of  the  distributed  magnetospheric  current  systems 
(Olson,  1974).  The  procedure  developed  for  obtaining  currents  is  quite 
different  from  previous  ones.  Instead  of  calculating  self-consistent  particle 
and  field  distributions  it  was  demanded  that  the  currents  produce  magnetic 
fields  similar  to  those  observed.  This  procedure  was  used  because  it  most 
readily  permits  the  development  of  quantitative  models  of  the  currents  and 
the  fields.  The  currents  that  produce  the  magnetic  field  features  described 
above  are  themselves  produced  by  the  drift  of  electrons  and  ions  in  opposite 
directions  in  the  non-uniform  total  magnetospheric  field  and  by  diamagnetic 
currents  caused  by  gradients  in  the  plasma  density.  In  the  inner  magnetosphere 
the  analysis  of  0G0-3  and  5 magnetometer  data  by  Sugiura  and  his  coworkers 
(Sugiura  et  al.,  1971;  Sugiura,  1972;  Sugiura  and  Poros,  1973)  provided 
important  information  on  the  currents  produced  by  trapped  particles.  In 
particular,  they  developed  contours  of  AB,  the  scalar  difference  between  the 
magnitudes  of  the  total  field  and  the  earth's  main  field.  Their  plots  provide 
an  indication  of  the  scalar  magnitude  of  the  field  produced  by  currents 
flowing  in  the  magnetosphere.  A noon  midnight  plot  of  AB  for  low  Kp  values  is 
shown  in  Figure  3-1.  The  considerable  structuring  in  the  equatorial  region 
suggests  that  currents  must  flow  there. 

The  procedure  for  determining  the  currents  was  then  as  follows.  The  general 
form  of  the  currents  as  suggested  by  the  observations  were  used  to  select  the 
positions  of  “wires"  and  the  strength  of  the  current  flowing  through  each  of 
them.  Integration  over  these  currents  then  yielded  values  of  the  field  from 
the  distributed  currents  that  were  added  to  the  magnetopause  magnetic  field 
(which  had  been  previously  determined)  and  the  earth's  main  field  and  com- 
pared with  observations.  Differences  in  the  observed  and  computed  values  of 
the  field  were  noted  and  changes  made  in  the  positions  and  strengths  of  the 
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currents  in  order  to  decrease  the  differences.  The  currents  in  the  inner 
magnetosphere  and  in  the  tail  were  considered  separately.  The  AB  contours 
obtained  from  the  final  model  are  shown  in  the  noon-midnight  meridian  plane 
in  Figure  3-2.  It  is  seen  that  the  model  correctly  reproduces  the  main 
observed  field  (older  models  produced  grossly  incorrect  AB  values). 

In  the  tail  regions  it  was  demanded  that  the  model  fit  both  the  observed 
decay  of  the  lobe  field  with  distance  down  the  tail  and  the  decay  of  the 
northward  component  through  th*>  equatorial  region  with  distance  down  the 
tail.  Several  attempts  were  made  at  fitting  these  features  with  the 
standard  “0"  current  loop  (e.g.,  Axford  et  al.,  1965;  Siscoe,  19f,6).  Witn 
the  cross  tail  currents  constrained  to  a thin  equatorial  sheet  it  was  not 
possible  to  model  both  the  lobe  and  equatorial  field  simultaneously.  A good 
fit  was  obtained  only  when  the  cross  tail  currents  were  allowed  to  flow 
through  the  plasma  sheet  and  the  return  currents  were  not  constrained  to  flow 
in  planes  of  constant  XSM.  In  order  to  fit  these  features  it  was  also 
necessary  to  have  the  loop  current  contribute  locally  to  the  northward  com- 
ponent of  the  field  at  the  equator.  The  observed  flare  of  the  radius  of  the 
tail  was  also  input  to  the  study.  Although  tests  with  varying  flare  rates 
indicated  that  the  tail  shape  had  only  a small  effect  on  the  field,  an 
analytic  expression  for  the  tail  flare  was  developed  from  reported  values 
(Siscoe,  1971)  for  use  in  the  study. 

Although  the  main  purpose  of  this  work  on  distributed  currents  was  to  provide 
an  input  for  the  development  of  a quantitative  magnetic  field  model,  it  also 
shed  light  on  several  other  magnetospherlc  features.  (1)  It  became  evident 
in  the  study  that  dynamic  processes  must  be  involved  in  the  maintenance  of 
the  quiet  magnetosphere.  The  magnetosheath  plasma  continuously  supplies  the 
distributed  currents  that  flow  across  the  tail  but  the  quiet  time  ring  cur- 
rent is  formed  by  trapped  particles  whose  drift  shells  do  not  intersect  the 
magnetopause.  Thus,  the  magnetosheath  cannot  directly  be  the  source  of  the 
inner  loop  currents.  They  must  instead  be  fed  by  the  tail  at  irregular 
(substora  or  storm)  intervals  and  so  can  be  formed  and  maintained  only  by 
temporal  changes  in  the  tail  magnetic  field  or  the  presence  of  a cross  tail 
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electric  field,  either  of  which  can  act  to  accelerate  particles  into  orbits 
that  close  around  the  dayside  magnetosphere.  (2)  A study  of  the  distant 
tail  field  suggests  that  the  currents  flowing  across  the  tail  may  be  larger 
near  the  boundary  of  the  plasma  sheet  but  are  not  particularly  strong  in  the 
neutral  sheet  region  and  that  the  neutral  sheet  should  be  thought  of  merely 
as  that  region  where  the  total  magnetic  field  is  quite  small.  (3)  These 
distributed  currents  also  contribute  to  the  magnetic  field  at  the  earth's 
surface,  and  therefore  must  be  included  in  the  study  of  the  quiet  daily 
variation  In  the  surface  field  (Sq)  and  in  the  study  of  D$t  which  is  a measure 
of  the  strength  of  the  ring  current.  These  magnetospherlc  currents  make  a 
significant  contribution  to  Sq  (the  observed  noon-midnight  difference  in  the 
Sq  field  at  mid-latitudes  Is  about  20  to  40  gammas  during  quite  times).  At 
the  equator,  the  noon-midnight  difference  in  the  north  component  of  the  field 
from  these  currents  is  just  under  9 gaimias  and  about  14  gammas  when  the 
effects  of  Induced  currents  are  Included.  However,  it  Is  still  believed  that 
they  are  not  the  predominant  source  of  Sq  as  was  suggested  by  Sarabhal  and 
Nalr  (1969)  and  that  Ionospheric  currents  produce  most  of  the  observed 
pattern  although  the  magnetospherlc  currents  must  contribute  significantly  to 
the  dally  variability  in  $q.  (4)  Magnetosheath  plasma  flows  across  the  tail 

of  the  magnetosphere  and  contributes  simultaneously  to  the  plasma  sheet  and 
the  distributed  currents  In  the  tall.  This  penetration  of  magnetosheath  plasma 
suggests  that  the  magnetosphere  Is  "open"  to  low  energy  particles  over  several 
regions  of  the  magnetopause.  This  work  suggested  that  both  high  and  low 
particles  can  enter  the  magnetosphere  even  If  It  Is  magnetically  closed.  Thus, 
although  details  of  merging  geomagnetic  and  Interplanetary  fields  have  not 
been  studied  quantitatively,  It  Is  quite  certain  that  the  magnetosphere  Is 
always  open  to  the  entry  of  charged  particles  and  that  In  any  self-consistent 
quantitative  model  of  magnetospherlc  particles  and  fields,  the  fields  will 
not  be  completely  confined. 

The  total  magnetospherlc  magnetic  field  can  be  found  by  integrating  over 
these  dlstriouted  currents  and  adding  the  field  value  to  the  field  produced 
by  the  magnetopause  currents  and  the  earth's  main  field.  However,  this  is 
a very  time  consuming  process.  For  the  model  to  be  useful  It  must  be  capable 
of  returning  the  total  magnetospherlc  magnetic  field  at  any  position  In  the 


magnetosphere  with  a minimal  expenditure  of  computer  time.  Therefore,  the 
values  of  the  magnetic  field  were  found  at  many  points  by  direct  integration 
over  the  distributed  and  magnetopause  currents.  These  values  of  the  field 
were  then  input  to  a least  squares  best  fit  computer  program  to  determine 
coefficients  of  the  polynominal  series  to  be  used  in  the  analytic  field 
representation.  Previous  representations  of  the  magnetopause  magnetic  field 
have  been  given  in  the  terms  of  the  scalar  potential.  Since  the  field  from 
the  distributed  currents  must  be  fit  in  the  region  of  the  source  currents, 
a scalar  potential  can  no  longer  be  used  since  V x B must  be  permitted  to  be 
non-zero.  Another  constraint  on  the  form  of  the  series  representation  is 
dictated  by  the  rather  sharp  structuring  that  the  field  from  the  distributed 
currents  exhibits  in  the  inner  magnetosphere.  Thus,  the  series  must  simul- 
taneously fit  a very  structured  region  near  the  earth  and  the  weak  field 
region  in  the  distant  tail.  A standard  one-dimensional  orthonormal  least 
squares  fitting  program  was  generalized  for  use  to  find  coefficients  in  n 
variables  (K.  A.  Pfitzer,  unpublished  reference,  1973).  In  order  to  fit 
both  the  structure  of  the  inner  magnetosphere  and  the  tail  field  topology, 
it  was  necessary  to  include  terms  of  up  to  sixth  order  in  a power  series 
plus  similar  terms  multiplied  by  an  exponential.  A final  requirement  for 
the  model  is  that  the  field  must  vary  smoothly  between  the  points  where  it 
was  computed  for  input  to  the  fitting  program.  Calculations  of  both  the 
distributed  magnetic  field  and  the  magnetopause  magnetic  field  show  that  the 
model  fields  do  vary  smoothly  over  the  region  where  they  are  defined.  This 
finding  is  particularly  important  for  the  distributed  maqnetic  field  since 
"wires"  were  used  initially  to  represent  source  currents  (Olson,  1974).  The 
form  of  the  series  representation  for  the  distribute*  lapnetic  field  and  the 
coefficients  have  been  published  in  Olson  and  Pfitzer  (1974).  A listing 
of  the  computer  deck  is  included  in  Appendix  A. 

Field  lines  at  various  latitudes  are  shown  for  different  local  times  in 
Figure  3-i  using  the  analytic  version  of  the  model.  Note  that  the  field 
lines  from  76°  magnetic  latitude  cross  the  equator  at  midnight  farther  from 
the  earth  than  the  82°  lines  shown  by  Head  (1964).  It  is  this  feature,  the 
extension  of  the  near  earth  field  lines  produced  by  the  inclusion  of  the 
distributed  magnetic  field,  that  makes  this  model  capable  of  reproducing 


many  particle  and  field  observations  that  previously  were  not  properly  modeled. 
Generally,  the  field  line  geometry  in  the  tail  remains  quite  dipolar.  There 
are  no  lines  that  run  parallel  to  the  equator  for  large  distances.  It  is 
emphasized  that  this  model  magnetic  fie’d  is  intended  to  represent  only  the 
quiet  (and  most  dipolar)  state  of  the  magnetosphere.  The  accuracy  of  the 
model  has  been  tested  by  comparing  model  calculations  with  particle  and  field 
observations.  Optical  tracking  data  (Adamson  et  al.,  1973)  from  the  NASA-MPE 
Barium  Cloud  Experiment  has  shown  that  the  field  lines  in  the  inner  magneto- 
sphere are  more  elongated  than  the  lines  calculated  from  models  consisting  of 
main  field  and  earlier  external  field  representations  (Barrish  and  Roederer, 
1973).  A computer  program  that  combines  any  series  expansion  of  the  main 
field  with  various  representations  of  the  magnetospheric  field  (K.  A.  Pfitzer, 
unpublished  manuscript,  1973)  used  the  GSFC  9/65  main  field  exDansion  and  the 
present  model  to  generate  field  lines  and  compare  them  with  the  optical  track 
of  the  barium  clouds  (see  Figure  3-4).  The  model  field  fits  the  observations 
quite  accurately  because  it  includes  the  depressed  inner  magnetosphere  field 
feature  produced  by  the  quiet  time  ring  current.  The  earth  intercept  of 
field  lines  from  geosnychronous  orbit  have  been  determined  for  various  local 
times  (see  Figure  3-5).  The  foot  of  the  field  line  lies  between  65  and  66° 
magnetic  latitude  depending  on  the  time  of  day.  This  is  almost  3°  lower  than 
the  values  given  by  previous  boundary  models  (Mead,  1964;  Olson,  1970).  The 
model  was  used  to  compute  high  and  low  energy  charged  particle  behavior  in 
the  magnetosphere.  The  high  latitude  earth  intercept  of  the  trapping  boundary 
calculated  from  the  model  agrees  well  with  the  observed  boundary  which  is 
four  to  five  degrees  lower  than  that  calculated  from  previous  models  (Pfitzer, 
1972).  The  trapping  boundary  computed  with  the  present  model  is  in  good 
agreement  with  observations  at  all  latitudes  (see  Figure  3-6).  The  high 
latitude  cutoff  for  cosmic  ray  particles  calculated  from  previous  models  have 
been  in  error  with  observation  by  5 to  7°.  With  the  present  models  the 
calculated  cutoffs  are  essentially  in  agreement  with  the  observed  values 
(Masley  et  al.,  1973).  It  is  the  depressed  field  region  from  3 to  6 radii 
geocentric  distance  produced  by  the  distributed  currents  that  causes  the  field 
lines  at  a given  latitude  at  the  earth's  surface  to  cross  the  equator  at 
larger  geocentric  distances  than  those  calculated  from  previous  models.  In 
turn,  it  is  this  field  geometry  that  allows  charged  particles  to  penetrate  to 
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lower  latitudes  than  is  possible  with  models  that  do  not  include  the  effects 
of  the  distributed  currents  that  flow  through  the  inner  magnetosphere. 

This  total  magnetic  field  model  has  been  restricted  in  two  ways.  First,  it 
represents  only  the  symmetric  magnetosphere.  This  geometry  occurs  only  when 
the  solar  wind  is  incident  perpendicular  to  the  earth's  magnetic  dipole  axis. 
The  second  "limitation"  of  the  model  is  that  it  represents  quiet  magnetic 
conditions.  Extensions  and  improvements  to  this  model  are  discussed  in 
Section  6.0.  (Work  on  the  development  of  this  field  model  was  supported  jointly 
by  this  contract,  the  Office  of  Naval  Research  (ONR),  and  the  MDAC  Independent 
Research  and  Development  Program.) 
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4.0  A MODEL  OF  THE  NEUTRAL  ATMOSPHERE 

Recent  satellite  experiments  and  previous  analyses  of  satellite  drag  data 
have  shown  that  the  high  latitude  neutral  atmosphere  Is  quite  structured. 

Both  the  Incoming  cusp  particle  fluxes  and  their  Influence  on  the  neutral 
atmosphere  have  been  measured.  To  understand  and  ultimately  predict  the 
effects  of  these  environmental  features  on  space  hardware  systems.  It  Is 
necessary  to  develop  a quantitative  model  of  the  density  of  the  upper 
atmosphere  that  Includes  the  energy  Input  from  the  dayside  cusp  particles. 

Such  a model  was  developed  under  this  contract  and  Is  described  in  this 
section.  It  Includes  both  the  solar  ultraviolet  and  charged  particle  energy 
sources.  The  model  Is  semi-empirical— It  Is  based  on  available  satellite 
data  and  our  present  understanding  of  both  the  UV  and  corpuscular  energy 
sources.  It  offers  a global  description  of  the  atmosphere  above  120  km. 

The  model  takes  Into  account  the  dependence  of  atmospheric  density  on  the 
value  of  the  solar  flux  constant.  The  UV  and  corpuscular  effects  are  com- 
puted In  separate  coordinate  systems  and  their  contributions  added.  It  Is 
believed  that  the  work  on  the  atmospheric  density  described  here  represents 
significant  progress  In  environmental  modeling.  It  Is  the  first  time  that 
corpuscular  energy  Inputs  have  been  considered  quantitatively  in  a neutral 
density  model.  It  Is  also  believed  that  this  Is  the  first  quantitative, 
analytic,  global  model  of  the  upper  atmospheric  density.  This  model  should 
be  useful  In  several  theoretical  studies  and  In  systems  applications. 

Ever  since  an  accel erometer  was  first  flown  on  a polar  satellite  (Bruce, 

1968),  It  has  been  apparent  that  there  are  at  least  two  density  bulges  In 
the  lower  thermosphere  during  geomagnetlcally  quiet  times.  The  low-latitude 
bulge  Is  usually  attributed  to  heating  by  solar  UV  radiation  (although  other 
energy  sources  may  contribute  Importantly),  but  the  high-latitude  bulge  Is 
almost  certainly  caused  by  energy  carried  by  the  plasma  which  makes  up  the 
solar  wind. 

There  are  several  charged  particle  sources  that  heat  the  upper  atmosphere,  but 
only  those  particles  that  precipitate  through  the  dayside  cusps  constitute  a 
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steady  source  with  corresponding  "permanent"  atmospheric  density  features. 
(Charged  particles  that  precipitate  Into  the  upper  atmosphere  during  sub- 
storms  and  magnetic  storms  produce  large  and  sudden  changes  In  density  but 
their  presence  Is  Irregular.)  The  dayslde  cusp  particles  consist  of  magneto- 
sheath protons  and  electrons.  The  protons  can  penetrate  to  about  110  km  and 
the  electrons  to  just  below  100  Ion  (Olson,  1971;  Olson,  1972).  The  protons 
carry  the  bulk  of  the  energy  that  Is  transferred  to  the  neutral  atmosphere. 
Their  peak  energy  loss  per  unit  change  in  altitude  is  at  169  km.  The  particle 
energy  transferred  to  the  atmosphere  eventually  produces  heating  which  causes 
neutral  particles  at  a given  altitude  to  rise  and  Increase  the  density  at 
heights  above  the  heating  source  region.  Thus,  It  is  expected  that  although 
atmospheric  density  should  be  Influenced  by  the  particles  to  altitudes  as 
low  as  100  km  the  largest  effects  should  occur  near  and  above  169  km. 

Two  coordinate  systems  are  used  In  the  r,.odel.  The  solar  UV  source  Is  best 
described  in  geographic  coordinates.  Thus,  Inputs  to  the  computer  subroutine 
must  include  the  universal  time  and  the  time  of  year.  Local  time  Is  also 
entered  as  a function  of  the  coordinates  of  the  point  where  the  atmospheric 
density  Is  to  be  determined.  The  dayslde  cusp  particles  are  constrained  to 
precipitate  Into  the  atmosphere  along  magnetic  field  lines.  The  cusp  Inter- 
section with  the  atmosphere  Is  nominally  15  degrees  below  the  magnetic  dipole 
axis  with  Its  longitude  center  on  the  magnetic  noon  meridian  plane  as  defined 
by  the  dipole  axis  and  the  sun-earth  line).  The  extent  of  this  Intersection 
In  longitude  is  about  12  hours.  The  extent  of  the  cusp  In  latitude  Is 
several  degrees  (3-5).  The  region  of  the  atmosphere  actually  heated  by  these 
particles  is  of  course  much  larger  because  the  Impact  energy  Is  spread  out  by 
winds,  thermal  conduction,  and  possibly  by  gravity  waves.  In  the  development 
of  the  present  model,  the  Input  energy  sources  have  been  considered  quantita- 
tively. The  determination  of  the  extents  of  their  heating  Influence,  however, 
has  not  considered  dissipation  mechanisms  such  as  winds  and  wave  phenomena. 
Rather,  the  extents  of  the  heating  regions  have  been  determined  from  satellite 
measurements  (both  direct  and  Indirect)  of  the  atmospheric  density.  The 
corpuscular  contribution  Is  then  found  in  magnetic  coordinates  and  combined 
with  the  UV  contribution.  The  form  of  the  equations  used  In  the  model  Is 
now  discussed. 


4.1  Model  Description 

In  the  Interim  Reports  to  AFOSR,  mathematical  descriptions  of  a neutral  density 
model  were  developed.  Those  reports  concerned  themselves  with  the  description 
of  the  functions  necessary  to  describe  the  neutral  density  model  and  the  devel- 
opment of  an  understanding  of  the  high  latitude  corpuscular  energy  source. 

Here,  a fast  FORTRAN  computer  deck  which  can  be  used  to  describe  the  atmos- 
phere from  0-450  km  Is  described.  In  several  cases,  shorter  expansions  or 
simpler  functions  than  discussed  In  the  Interim  Reports  were  used  whenever 
such  usage  did  not  compromise  the  accuracy  of  the  model  within  current  experi- 
mental error.  We  believe  that  during  quiet  times  In  the  range  120-450  km  (the 
principle  altitude  range  covered  by  this  study)  the  model  has  an  accuracy  of 
about  10  percent.  From  0-120  km  a simple  power  series  expansion  extends  the 
120-450  km  model  to  the  ground  by  fitting  It  to  the  CIRA  72  mean  reference 
atmosphere  (COSPAR,  1972).  The  average  deviation  from  CIRA  72  between  0-120 
km  Is  less  than  one  percent  with  a 7th  order  fit. 

The  FORTRAN  deck  has  as  inputs  position  and  time  as  well  as  parameters  des- 
cribing solar  activity.  The  input  parameter  for  the  UV  heating  Is  the  10.7  cm 
solar  flux  and  the  Input  parameter  for  the  particle  heating  Is  proportional 
to  the  particle  flux  and  the  location  and  size  of  the  cusps. 

The  following  paragraphs  define  the  equations  utilized  In  the  development  of 
the  FORTRAN  code.  The  entire  functional  representation  Is  contained  In  this 
report,  so  that  the  reader  can  develop  his  own  FORTRAN  code  if  so  desired. 
Appendix  B contains  the  MOAC  FORTRAN  code  for  the  density  model.  It  has 
been  carefully  optimized  to  minimize  the  calls  to  elementary  functions  (sines, 
cosines,  square  roots,  etc.)  and  to  minimize  the  number  of  multiplies.  A test 
program  and  a sample  computer  output  are  also  listed  to  assist  the  reader  In 
usage  of  the  program. 

4.1.1  Functional  Form  for  the  Total  Density 

The  total  density  Is  expressed  as  the  sum  of  two  terms:  the  first  term,  pp, 
describes  the  combined  effects  of  ultraviolet  heating  and  effects  other  than 


those  due  to  the  corpuscular  energy  sources.  The  second  term,  Apc,  gives 
the  contribution  due  to  the  density  associated  with  corpuscular  heating 
caused  by  particles  entering  the  atmosphere  through  the  magnetospherlc 
cusps.  Thus  we  may  write  the  density  as 

p ■ p(z,  X,  <j>,  T,  D,  F,  P)  B Py  + Apc  (4-1) 

where 

3 

pis  the  total  neutral  density  In  gm/cm 
Pp  Is  the  density  due  to  non-corpuscular  energy  sources 
Apc  Is  the  density  due  to  corpuscular  energy  sources 
z Is  the  altitude  In  km 

X Is  the  geographic  latitude  in  degrees  (+  Is  North,  - Is  South) 

<p  Is  the  geographic  longitude  In  degrees  east  (0-360°) 

T Is  the  universal  time  In  hours 
D Is  the  day  of  year 
F Is  the  declmetrlc  solar  flux 
P represents  the  particle  heating  parameters 

The  function  py  Is  expressed  as  a product  of  four  terms  (See  Eq.  (2)  Interim 
Report  Sept  75) 

Py  ■ P0(Z,  F)  B(z,  X,  t)  J(z,  X,  D)  Q(z,  D)  (4-2) 

where  t Is  the  local  time  In  hours  and  Is  given  by  t « T + <J>/15, 
p0  Is  the  mean  equatorial  density  as  a function  of  height, 

B gives  the  diurnal  variation  ,.\  the  density, 

J gives  the  seasonal  and  latitudinal  dependence, 

Q Is  the  semi-annual  variation. 

Each  of  the  above  terms  Is  now  explained  In  some  detail. 

4.1.2  Altitude  Dependence 

For  average  conditions  the  mean  equatorial  neutral  density,  p0(z,  F),  can  be 
expressed  as 


p0  * p( 1 20)  exp 


(4-3) 
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1 .737 /z  - 103 
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for  120  km  < z < 450  km, 


-11  3 

where  p(120),  the  mean  equatorial  density  at  120  km.  Is  2.7  x 10  gm/cm  . 
The  square  root  relationship  In  the  denominator  of  the  exponent  Is  used 
because  of  experimental  evidence  that  scale  heights  In  the  thermosphere 
vary  approximately  as  the  square  root  of  the  altitude  (Moe,  et  al.,  1968). 


At  the  lowest  altitudes,  p0  is  based  on  Moe  (1973).  At  the  higher  altitudes, 
parameters  are  chosen  to  make  this  model  In  approximate  agreement  with  the 
forthcoming  U.S.  Standard  Atmosphere. 


Equation  4-3  can  be  modified  to  take  Into  account  the  dependence  of  the 
equatorial  density  on  solar  activity  by  replacing  the  constant  1.737  with 
the  function  A,  where 

A * 0.99  + .518  (t^- 

F Is  the  10  cm  solar  flux  constant  measured  at  Ottowa,  Canada  and  F Is  Its 

m09 

average  over  the  preceding  three  months.  They  are  given  In  units  of  10 
watts/m2/Hz.  The  function  A has  been  developed  to  represent  the  solar  cycle 
variations  deduced  from  satellite  measurements  by  Klng-Hele  and  Quinn  (1966), 
and  Moe  (1969). 


1/2 


(4-4) 


4.1.3  Diurnal  Variation 

The  function  B (z,  X,  t)  which  represent  the  diurnal  variation  Is  given  by 
the  expression  B (z,  X,  t)  ■ [1  ♦ (f(t)  - 1)  cos  X]y^*  where  f(t)  Is  the 
diurnal  variation  In  density  at  the  equator.  The  function  f(t)  Is  repre- 
sented as  a Fourier  series  such  that  given  u(z)  ■ 1 at  450  km,  the  density 
at  the  equator  accurately  models  the  density  variation  observed  by  Hedln 
et  al.  (1974).  Development  of  the  series  Is  described  In  Appendix  D. 
However,  use  of  the  full  six  term  Fourier  series  expansion  Is  not  justified 
at  this  time  because  of  the  large  scatter  In  the  data  of  Hedln  et  al.  In 


30 


]-*  lWlfc 


■y* f V a 


the  FORTRAN  deck  we  limit  the  expansion  to  the  first  three  terms.  This 
modifies  the  fit  by  only  two  percent  and  Is  well  within  the  ten  percent 
accuracy  of  Hedln  et  al.'s  curve.  Thus 

f(t)  * .994  + .545  cosC^  - 14.745)]  + .102  cos  [J(t  - 1.838)]  (4-5) 

The  function  ^)  which  gives  the  height  dependence  of  the  diurnal  variation 
must  satisfy  the  following  conditions: 


1.  p(450)  5 1.  I.e.,  f(t)  must  reproduce  Hedln  et  al.'s  values 

at  450  km 


2.  p(186)  - .4  In  order  to  match  the  LOGACS  data  near  noon 


3.  p(230)  3 .57  In  order  to  fit  the  data  of  Chlng  and  Carter  (1974) 

near  local  midnight 

4.  y(120)  ■ 0.  I.e.,  the  diurnal  variation  must  approach  zero  at 

120  km,  the  lower  limit  of  the  model  where  It  Is  to 
Interface  with  CIRA  72, 


The  expression  y(z)  ■ 1.1  x (1  - exp(-(z-120)/150))  satisfies  the  above  four 
conditions. 


4.1.4  Latitudinal  and  Seasonal  Dependence 

The  1974  and  1975  Interim  Reports  differ  substantially  In  their  treatment  of 
the  term  J(z,  X,  0)  which  describes  the  latitudinal  and  seasonal  dependence. 
The  1975  Interim  Report  Introduces  a relatively  complex  function  which  depends 
on  18  empirically  determined  coefficients  In  order  to  Improve  the  fit  of  the 
density  model  with  the  mid-latitude  LOGACS  data.  Although  this  function 
significantly  Improved  the  fit  to  the  LOGACS  data.  It  Introduced  a rather 
large  seasonal  latitude  dependence.  For  example,  at  solstice  a ratio  In 
the  North  to  South  polar  densities  of  10  to  1 Is  Introduced  at  186  km.  At 
450  km  the  ratio  exceeds  100  to  1.  Hedln  et  al.  (1974)  observe  only  a 2 to  1 
latitudinal  variation  at  450  km.  Therefore,  even  though  the  more  complex 
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function  does  improve  the  fit  to  LOGACS  during  northern  summer  day,  it  cannot 
be  used  at  this  time  in  a FORTRAN  deck  representing  a global  density  model. 
Therefore,  the  much  simpler  expression  defined  in  the  1974  Interim  Report  is 
used  in  the  model  construction.  We  thus  write 

my(z) 

J(z,  X,  D)  = [1  - a(l  - cos(X  - X$))]  (4-6) 

As  is  the  latitude  of  the  sun's  subsolar  point  and  is  given  by  As  = 

23.5  sin  — • A value  of  a = .25  is  consistent  with  observations 
indicating  only  a small  latitude  dependence  due  to  UV  heating.  A value  of 

2.5  for  m Is  also  consistent  with  the  requirement  that  rru(186)  a 1 at  186  km 
to  best  fit  the  LOGACS  data  and  mp(450)  = 2.5  in  order  to  best  fit  the  lati- 
tude dependence  observed  by  Hedin  et  al.  (1974)  at  450  km.  Since  mu(z) 
approaches  zero  as  z approaches  120  km,  the  function  J approaches  unity  as 

z approaches  120  km,  and  thus  the  latitude  and  seasonal  variation  goes  to 
zero  as  z approaches  120  km. 

4.1.5  The  Semi-annual  Variation 

Derivation  of  the  semi-annual  variation  represented  as  the  term  Q In  Eq  (4-2) 
Is  explained  In  considerable  detail  In  Appendix  E.  The  function  Q is 
expressed  as 

Q(z,  D)  « 1 + [R(z)  - 1]  G(D)  (4-7) 


where 


R(z)  - .98  + .27  x 10’2z  - .85  x 10"6  z2  - .59  x 10"9  z3  (4-8) 


and 


G(D)  « .143  cos(§jk<D  - 4))  + .239  cos(^<D  - 109))  (4-9) 

+ .044  cos(|Jj(D  - 66)) 
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The  semi-annual  variation  does  not  approach  zero  as  z approaches  120  km. 
Observations  Indicate  that  a significant  diurnal  variation  Is  observed  only 
to  as  low  as  90  km.  Therefore,  Q Is  allowed  to  vary  below  120  km  such  that 
no  diurnal  variation  exists  below  75  km  and  the  total  density  agrees  with 
CIRA  72.  This  damping  is  explained  In  Section  4,1.7  which  described  the 
model  from  0-120  km. 

4.1.6  Corpuscular  Heating 

Apc,  the  density  Increase  due  to  the  corpuscular  heating,  is  given  In  the  1974 
Interim  Report  as 

Ape  = Py  * C(z,  Am,  tfl,  P)  (4-10) 

and  In  the  1975  Interim  Report  as 

Ape  * Po  C(z,  Ajjj,  tfij,  P)  (4-11) 

where  C Is  the  cusp  heating  function.  We  have  chosen  to  use  equation  4-10 
Instead  of  4-11  because  studies  by  Hedln  et  al..  Indicate  a definite  high 
latitude  peak  In  the  summer  hemisphere,  but  none  or  a smaller  one  In  the 
winter  hemisphere.  Equation  4-10  Introduces  a summer-winter  dependence 
Into  the  cusp  heating  term  whereas  equation  4-11  has  north  south  particle 
heating  symmetry  for  all  seasons.  The  term  C Is  given  a scale  height  factor 
of  about  50  km  since  we  wish  to  have  zero  heating  at  120  km  and  approach 
maximum  heating  near  170  km  (Olson;  1971,  1972).  Thus 

C(z,  Am,  tm,  P)  - exp(-  (z  - 120)/50)C(Am,  t*.  P)  (4-12) 

The  function  C(Am,  t^,  P)  Is  dependent  on  the  magnetic  latitude.  An,,  the 
magnetic  local  time,  tm,  and  the  cusp  heating  parameters  represented  In  the 
above  equation  by  P.  The  parameters  represented  by  P are  the  particle  pre- 
clpatatlon  Intensity,  I,  the  half  width  of  heating  region  In  magnetic  latitude, 
e0,  the  half-width  of  the  heating  region  In  magnetic  local  time,  E,  and  the 
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magnetic  latitude  of  the  center  of  the  northern  cusp  boating  region,  Amc. 

The  geomagnetic  coordinate  system  Is  defined  such  that  the  North  geomagnetic 
pole  is  located  at  69°  W and  78.5°  N geographic.  The  zero  in  magnetic 
longitude  is  along  69°  W geographic  longitude. 


Magnetic  local  time,  tm,  Is  defined  such  that  zero  is  in  the  anti -solar 
direction.  The  magnetic  local  time,  tm*  Is  measured  in  degrees.  The  magnetic 
local  time  of  the  observation,  tm0.  is  given  by 


tmo  " 15  * UT  + <j>m  - 69 


(4-13) 


where  is  the  magnetic  longitude.  The  local  time  of  the  center  of  the  cusp, 
tmc,  is  always  located  at  the  subsolar  longitude  (magnetic  local  noon)  and 
thus,  tmc  * 180°.  The  difference  In  magnetic  local  time  between  the  observa- 
tion point  and  the  cusp  center,  Atm,  Is  given  by 


Atm  ■ 15  UT  + $m  - 69  - 180 


- 15  UT  + <fcn  - 249 


A^  is  restricted  to  the  range  -180°  to  +180°.  Therefore 


(4-14) 


if  Atm  > 180 
If  Atm  < 180 


Atm  - - 360 

Atm  - Atm  ♦ 360 


The  magnetic  latitude,  Am.  is  restricted  to  i 90°  (+  Is  North  - Is  South). 

The  difference  in  magnetic  latitude,  AXm,  between  the  observation  point,  Amo* 
and  the  center  of  the  cusp,  AmCt  Is  given  by  the  equation 


AXm  ■ |Xmol  * IXfflcI 


(4-15) 


Given  the  above  definitions  we  can  now  express  the  function  C by  the 
following  expressions: 
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C<V  V p>  ■ 0 
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C(Am,  t , P)  ■ 1 x L * (1  + cos  IjE- ) 


If  l«J  < s 


(4-16) 


where 


At. 


L = 2 - 


m' 


Ae  * AA 


If  Ltrr  i * 
m 


(4-17) 


L * 1 


Ae  - cos^sInA^sInA^+cosA^cosA^cos  (fAAj-E)] 


If  At  > E 
m 


The  above  defined  expression  for  C defines  a density  bulge  that  wraps  itself 
around  the  geomagnetic  pole  and  Is  most  Intense  at  local  magnetic  noon. 


4.1.7  Density  Below  120  km 

In  order  to  make  the  FORTRAN  deck  a more  useful  tool,  a fit  to  the  CIRA  72 
mean  reference  atmosphere  was  added  to  the  model.  The  junction  between  the 
0-120  km  model  and  the  120-450  km  model  Is  everywhere  continuous  and  the 
slopes  are  closely  matched,  but  not  mathematically  continuous.  Below  120  km 
the  density  Is  given  by 


p - (1.2252  x 10‘3  exp(  E7  a,z1))Q' 

1-1  1 


(4-18) 


where 


Q*  ■ Q when  z > 90  km 


Q'  ■ 1 ♦ (Q  - 1)  exp(-(z-90r/200)  when  z < 90  km 


(4-19) 


Q Is  the  diurnal  variation  defined  by  equation  4-7.  The  coefficients  for 
the  expansion  In  equation  4-18  are 
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a1  * -.8015593  x 10'1 
a2  « -.3366924  x 10"2 
a3  = +.1821500  x 10‘4 
a4  = +.1739634  x 10’5 
a5  = -.3658010  x 10‘7 
ag  = +.2688846  x 10'9 
a?  * -.6820457  x 10‘12 


The  seventh  order  expansion  produces  a fit  to  CIRA  72  below  75  km  which  has 
an  RMS  error  of  less  than  one  percent  and  a maximum  error  of  2.5  percent. 
Additional  terms  can  easily  be  added  to  Improve  the  fit. 


4.2  Model  Results  and  Comparisons  with  Observations 

This  model  of  the  density  of  the  neutral  upper  atmosphere  is  functionally 
very  simple  and  computationally  very  fast.  Yet,  It  considers  most  of  the 
known  variations  In  the  atmosphere.  It  treats  latitude,  longitude,  diurnal, 
seasonal,  semi-annual,  altitude,  solar  cycle  and  particle  precipitation 
variations  and  can  thus  be  used  to  predict  the  density  of  the  atmosphere  as 
a function  of  time,  position,  solar  and  magnetospherlc  parameters.  Since  It 
Includes  the  heating  effect  of  particles  precipitating  Into  the  high  latitude 
regions.  It  Is  especially  valuable  In  the  study  and  prediction  of  the  density 
of  the  polar  atmosphere.  Furthermore,  since  this  high  latitude  effect  Is 
parameterized.  It  can  be  used  to  predict  Increases  In  the  high  latitude 
density  when  the  solar  wind  particle  flux  Increases  and  the  cusp  location  Is 
observed  to  change.  The  model  Is  analytic  and  thus  differentiable.  Only 
three  regions  exist  where  the  derivative  of  the  density  Is  discontinuous. 

These  are  at:  120  km  where  the  120-450  km  model  is  matched  to  CIRA  72;  over 
the  geographic  poles;  and  over  the  geomagnetic  poles.  In  all  three  cases, 
however,  the  density  Is  continuous  and  the  slopes,  although  not  mathematically 
continuous,  are  closly  matched. 

The  three  dimensional  mercator  projections  In  Figures  4-1,  4-2  and  4-3  give 
an  overview  of  the  model's  features.  Figure  4-1  shows  the  density  at  400  km 
during  summer  solstice  at  a universal  time  of  1600  hours.  The  UV  density  peak 
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FIGURE  4-1.  Density  at  400  km  near  sumner 
solstice  at  a universal  time  of 
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FIGURE  4-2 
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Is  located  to  the  north  of  the  equator  and  at  about  1400  hours  local  time. 

The  particle  heating  peak  Is  most  Intense  on  the  subsolar  magnetic  longitude 
containing  the  magnetic  dipole.  Since  the  total  density  Is  dependent  on  the 
particle  heating  effect  as  well  as  the  UV  heating,  the  northern  particle 
heating  peak  which  Is  In  sunlight  is  much  more  pronounced  than  the  southern 
particle  heating  peak  which  Is  far  into  darkness  during  the  summer  solstice 
night.  Notice,  furthermore,  that  the  North  pole-South  pole  density  differ- 
ence during  solstice  Is  about  a factor  of  2 and  Is  In  agreement  with  Hedln 
et  al.  (1974).  Figure  4-2  Is  the  same  as  4-1  except  that  the  universal  time 
is  equal  to  zero.  The  southern  particle  heating  peak  seen  In  Figure  4-2  is 
now  somewhat  larger  since  the  southern  dipole  Is  now  on  the  day  side  of  the 
earth.  The  UV  heating  peak  still  lags  the  sun  by  about  two  hours.  Figure 
4-3  Is  a similar  mercator  projection  at  150  km.  At  150  km,  the  amplitudes 
of  the  UV  and  the  particle  heating  peaks  are  considerably  reduced.  All  of 
the  variations  except  the  semi-annual  variation  approach  zero  as  the  altitude 
approaches  120  km.  The  semi-annual  variation  approaches  0 as  the  altitude 
approaches  75  km. 

A three  dimensional  polar  plot  Is  used  In  Figure  4-4  to  represent  the  density 
bulge  introduced  at  high  latitudes  due  to  particle  precipitation.  The  particle 
heating  region  surrounds  the  geomagnetic  pole  and  Is  maximum  In  the  solar 
direction  and  drops  to  zero  In  the  anti-solar  direction.  Its  amplitude, 
location  In  geomagnetic  latitude,  and  Its  size  are  adjustable  parameters 
related  to  solar  wind  and  magnetospherlc  properties.  The  adjustable  param- 
eters have  been  set  to  give  a best  fit  to  the  L0GACS  data  during  quiet  times. 
Figure  4-5  shows  a comparison  between  the  current  model  and  the  L0GACS  data. 

The  current  model  Is  able  to  represent  both  the  low  latitude  UV  bulge  as  well 
as  the  high  latitude  particle  heating  bulge.  Several  passes  of  L0GACS  data 
were  available  and  these  were  used  to  determine  the  particle  heating  intensity 
parameter,  I.  I ■ .15  gives  a best  fit  to  the  L0GACS  data.  It  Is  the  param- 
eter I that  must  be  varied  to  compensate  for  changes  In  the  precipitating 
particle  fluxes.  The  static  diffusion  models  (e.g.,  Jacchla,  1971)  do  not 
model  this  high  latitude  phenomenon. 
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LATITUDE  (DEG) 

FIGURE  4-5.  Comparison  of  the  model  with  the 
data  of  the  LOGACS  satellite. 

The  LOGACS  data  has  been  reduced 
to  186  km.  F a F * 160  and 


The  semi-annual  variation  was  adjusted  to  match  the  results  of  Hedln  et  al. 
(1974)  at  450  km.  Figure  4-6  shows  the  comparison  between  the  current  model 
and  Hedln  et  al.'s  0G0-6  data.  The  fit  at  450  km  for  both  0 and  45  degrees 
Is  within  experimental  error.  The  amplitude  of  the  diurnal  variation 
decreases  with  decreasing  altitude.  The  rate  of  decrease  was  adjusted  to 
fit  the  LOGACS  data  (Figure  4-5)  at  186  km  near  local  noon  and  the  data  of 
Ching  and  Carter  (1974)  at  230  km  near  local  midnight  (Figure  4-7).  the 
data  of  Ching  and  Carter  were  taken  during  April  1972  when  the  10  cm  solar 

Op  p 

flux  was  115  x 10  watts/m  /Hz.  The  agreement  In  absolute  density  is 
excellent  except  for  high  northern  latitudes.  For  night  time  and  latitudes 
above  60  degrees,  this  model  as  well  as  others  consistently  underestimate 
the  density  by  about  20  percent.  This  region  on  April  24th  Is  very  near  the 
light-dark  terminator  and  perhaps  the  model's  dependence  on  the  rate  of  change 
of  density  with  respect  to  the  solar  angle  is  too  rapid.  The  effect  Is  pro- 
bably not  due  to  neglected  particle  precipitation  effects  since  similar 
discrepancies  do  not  occur  near  the  south  pole.  The  southern  measurements 
were  made  well  within  the  darkness  region  and  were  far  removed  from  the 
light-dark  terminator  over  the  pole.  This  region  where  the  difference 
between  observation  and  model  averages  about  20  percent  Is  small  In  size. 

In  almost  all  other  regions  where  the  model  has  been  tested,  agreement 
between  model  and  data  Is  better  than  10  percent. 

Several  additional  comparisons  between  the  model  and  average  atmospheric 
densities  can  be  performed  by  averaging  outputs  from  the  model.  Figure  4-8 
shows  the  altitude  dependence  of  the  density  for  solar  minimum  (F  ■ F » 70) 
and  solar  maximum  (F  = F ■ 220)  at  0400  and  1400  hours  local  time  (l.e., 
near  the  minimum  and  the  maximum  In  the  diurnal  cycle).  Also  shown  i-re  points 
from  the  CIRA  72  mean  reference  atmosphere  which  represent  a diurnal ly  and 

seasonally  averaged  atmosphere  at  a latitude  of  30  degrees  when  the  10  cm 

-22  2 

solar  flux  has  a value  of  145  x 10  watts/m  /Hz.  The  dashed  curve  repre- 
sents the  model  with  F ■ F * 145  and  averaged  over  universal  time  and  the 
days  of  the  year  at  a latitude  of  30  degrees  North  and  a longitude  of  0 
degrees  East.  We  note  the  excellent  agreement  with  CIRA  72. 
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FIGURE  4-7.  The  nighttime  data  of  Chlng  and 
Carter  (1974)  Is  compared  with 
the  model.  The  average  10  cm 
solar  flux  during  April  1972  was 
^ 115.  This  value  was  input  to 
the  model  for  the  comparison. 


FIGURE  4-8.  Density  versus  altitude  for  solar 
minimum  (F  ■ F * 70)  and  solar 
maximum  (F  ■ F ■ 220)  near  the  dally  maximum 
(IT  ■ 1400  hrs)  and  the  dally  minimum  (IT  ■ 

0400  hrs).  Also  shown  Is  a comparison  between 
a universal  time  and  the  seasonally  averaged 
model  output  at  a latitude  of  30*  (F  ■ F ■ 145) 
with  the  CIRA  72  mean  reference  atmosphere. 
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Hedln  et  al.  (1974)  also  presents  a diurnal ly  averaged  curva  for  solstice 
showing  the  data  as  a function  of  latitude.  When  the  present  model  is  aver- 
aged over  universal  times  and  geographic  longitudes  at  solstice  a curve  very 
similar  In  appearance  to  Hedln  et  al.'s  060-6  curve  is  obtained  (See  Figure 
4-9).  In  the  summer  hemisphere  a second  high  latitude  maximum  (separate  from 
the  UV  maximum)  is  produced  due  to  the  particle  precipitation.  This  agrees 
well  with  the  double  maximum  observed  by  060-6.  The  major  differences  between 
the  060-6  data  and  the  model,  defined  by  this  report,  are  that  the  060-6  data 
are  about  25  percent  higher  and  that  the  model  determined  summer  high  latitude 
peak  Is  slightly  low  (when  the  L06ACS  determined  value  for  I,  I = .15,  is 
used).  Hedln  et  at.  state  that  their  seasonal  dependence.  In  particular  the 
annual  and  semi-annual  variation  Is  poorly  defined  and  thus  the  absolute 
amplitude  has  a large  uncertainty.  The  shape  In  latitude  dependence  is, 
however,  quite  accurate.  The  sumner  hemisphere  high  latitude  difference  is 
easily  corrected  by  Increasing  the  particle  heating  parameter  by  - 30  percent 
to  I » .2.  Since  no  data  are  available  during  this  time  period  for  determining 
the  average  particle  flux  entering  the  cusp,  a value  of  .2  Is  not  Inconsistent. 
The  average  variation  predicted  by  the  model  Is  In  substantial  agreement  with 
the  060-6  data. 

4.3  Use  of  Model  to  Describe  the  Disturbed  Atmosphere 
The  model  described  above  was  developed  with  a very  general  set  of  Input 
parameters  that  describe  variability  In  charged  particle  and  ultraviolet  heat 
sources.  To  date  those  parameters  have  been  adjusted  such  that  the  model 
describes  the  "undisturbed"  atmosphere.  Work  on  the  description  and  parameter- 
ization of  the  disturbed  atmosphere  was  initiated  In  the  past  year  and  described 
in  the  1975  Interim  Report.  By  appropriately  determining  the  dependence  of  the 
input  parameters  on  magnetic  and  solar  activity  and  taking  into  account  storm 
time  dependent  behavior,  the  model  can  be  extended  for  use  over  a range  of 
disturbance  conditions.  For  example,  charged  particle  heating  during  magneto- 
spheric  substorms  (via  the  plasma  sheet  lobes)  and  during  magnetic  storms  (from 
the  decay  of  the  ring  current)  can  readily  be  Incorporated  into  the  model  by 
appropriate  expansion  of  the  dayslde  cusp  region.  (Particle  precipitation 
through  the  lobes  of  the  plasma  sheet  will  extend  the  corpuscular  heating  region 


to  all  longitudes  and  to  lower  altitudes  on  the  nightside  due  to  the  relatively 
larger  primary  particle  energies  - from  10  to  50  keV.  Particle  heating  from 
ring  current  decay  will  extend  the  corpuscular  energy  source  to  much  lower 
latitudes  and  introduce  a dependence  on  storm  onset  time.)  As  has  been  our 
general  modeling  philosophy  (see  Section  6.0)  the  variability  in  the  input 
parameters  could  be  determined  by  comparing  preliminary  model  output  data  with 
observations.  As  with  our  magnetic  field  modeling  work,  an  attempt  is  made  to 
represent  all  conditions  with  one  basic  model.  For  purposes  of  prediction  and 
satellite  monitoring  we  also  believe  that  the  input  parameters  should  be 
capable  of  direct  and  continuous  observation  and  not  be  temporally  and  spatially 
averaged  indicies. 


5.0  SOLAR  WIND  AND  NEAR  FARTH  ENVIRONMENTAL  JgAMgERS 

The  models  discussed  in  the  previous  two  sections  (end  all  oth" 
environmental  models)  require  input  parameters  that  can  be  monitored  ” ™ 
tin*  if  the  models  are  to  be  used  for  prediction  purposes.  os- 
made  of  several  observational  data  sets  to  find  observables  that  directly 
described  near  earth  environmental  features  (e.g..  the  extent  of  davsrde 
aurora  provides  a gross  picture  of  the  location  of  the  dayside  ™sps). 

Possible  correlations  between  solar  wind  parameters  and  near  earth  en  r " - 
mental  features  were  also  sought.  (They  would  provide  a source  P 
input  parameters  for  the  models.) 

Density  measurements  deduced  from  orbital  perturbations  of  the  balloon 
satellites  Explore  19  and  Explorer  24  (dacchia  and  Slowey,  1972)  were  granhed 
fur  comparison  with  Vela  3 and  5 solar  wind  measurements  over  a peno  o 
several  years  in  order  to  determine  how  the  relationships  vary  « » 
soot  cycle  Comparative  data  covering  four  months  are  shown  in  Figures 
and  Illustrate  the  relationships  for  a wide  range  of  solar  anMejaqnet  < 

activity  The  correlation  between  air  drag  and  solar  wind  soee 
Suggests  that  nearly  every  change  in  the  solar  wind  produced  a corres- 
ponding change  in  the  energy  input  to  the  upper  atmosphere. 

Evidence  for  a high-latitude  density  bulge  at  guiet  timer . was  furnished  bv  the 

low-g  accelerometer  aboard  the  Air  Force  satellite  v 

(Bruce , 1968.  1972;  Oe  Vries.  1972).  In  Figure  6*  selected  I GA  data 

replotted  to  show  two  energy  sources  not  explicit  y in.  u e 

diffusion  -dels.  The  solid  line  in  the  Figure  represents  the  o 

density  when  the  geomagnetic  planetary  amplitude.  Ap  was  , in 

conditions.  (The  nocturnal  values  of  density  are  extrapol.  edj  he  a 

were  adjusted  to  a constant  altitude  of  100  nautical  H«  86  k”> ' e„ 

observed  latitudinal  variation  of  density  differs  considerably  from  that 

b,  the  dacchia  Model  (1970).  The  measured  density  appears  tc ’ ^treviolet 
expansion  In  response  to  two  discrete  energy  sources:  0)  The 

source  at  low  latitudes  and  (2)  the  corpuscular  source  at  high. 
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5 . 1 Composition  Measurements 

Longitudinal  variations  of  thermospheric  composition  have  been  observed  bv 
the  mass  spectrometer  on  0G0  6 (Hedin  and  Reber,  1972),  who  attribute  them 
to  geomagnetic  control  of  polar  heat  input.  The  contours  of  concentration 
in  the  Antarctic  atmosphere  at  two  local  times  are  illustrated  in  Figures  5-6 
and  5-7.  Their  geographic  extent  correlates  well  with  the  average  cusp 
location.  The  absolute  maximum  Ng  concentration  occurs  at  about  0800  LIT 
(3  hours  after  local  noon)  at  an  invariant  latitude  of  70°S.  Similar  per- 
sistent longitudinal  variations  of  positive  ions  occur  (Taylor,  1971).  (It 
requires  some  time  for  the  density  and  composition  to  respond  to  energy 
sources,  while  a nearly  instantaneous  response  is  provided  by  the  airglow.) 

5.2  Airglow  and  Auroral  Measurements 

O 

The  ISIS- I I Satellite  generated  global  maps  of  the  6300  Angstrom  line  of  atomic 
oxygen  (Shepherd  and  Thlrkettle,  1973).  The  most  prominent  feature  observed  is 
a band  of  permanent  red  aurora  on  the  dayslde  of  the  earth,  centered  on  magnetic 
noon  at  about  78°  invariant  latitude.  The  intensity  contours  of  this  aurora 
measured  near  0535  UT  on  14  December  1971  are  shown  in  Figure  5-8.  The  over- 
all shape  resembles  the  contour  of  the  dayside  cusp  at  0535  UT.  The  computed 
instantaneous  cusp  region  is  projected  on  the  airglow  contours  in  Figure  5-9. 

The  airglow  measurements  furnish  detailed  information  about  the  instantaneous 
energy  distribution  within  the  cusp.  Figure  5-8 also  shows  the  aurora  on  the 
nlghtslde,  which  is  produced  by  precipitation  from  the  tail.  By  coidbining 
the  ISIS- I I airglow  observations  with  magnetospheric  modeling  calculations, 
information  about  the  atmospheric  energy  sources  provided  by  the  magnetosphere 
during  the  day  and  night  at  quiet  and  disturbed  times  can  be  obtained. 

Measurements  of  Lyman-alpha  aurora  and  airglow  have  been  reported  by  Clark  and 
Metzger  (1969)  and  by  Metzger  and  Clark  (1971),  who  flew  a narrow-band 
scanning  photometer  on  the  Air  Force  research  satellite,  0V1-10.  This  experi- 
ment clearly  showed  the  enhancements  which  occur  during  geomagnetic  storms, 
as  well  as  the  variation  with  latitude  during  quiet  times.  In  the  various 
measurements  reported,  the  inferred  particulate  flux  in  the  auroral  zone  varied 
from  less  than  1 erg/cm  sec  during  several  quiet  periods  to  36  at  one  point 
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FIGURE 


-6  C 0 II  T U U R PLOT  I II  GEOGRAPHICAL 

COORDINATES  OF  N 2 DENSITY 
EXTRAPOLATED  TO  AGO  K II  FOR 
2OO0-2AOO  LT  AND  AVERAGED 
OVER  DATA  PROM  AUGUST  23  TO 
SEPTEMbEK  A , 19G1  WHEN  A p — 1 2 . 
DASHED  LINES  INDICATE  LACK 
OF  DATA.  DP  INDICATES  DIP  POLE 
AND  G H INDICATES  GEOMAGNETIC 
POLE. 


(AFTER  II  E D I N AND  R E B E R ) 
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180* 


figure  5-7  CONTOUR  PLOT  OF  H 2 DENSITY 
UNDER  THE  SAME  CONDITIONS 
AS  FIGURE  3-7  EXCEPT  FOR 
0A00-080U  LT. 


(AFTER  IIEDIN  AND  REBER) 
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figure  5-9  COMPUTED  REGION  OF  PRECIPI- 

TATION FROM  THE  D A Y S I D E 
CUSP  (SHADED)  COMPARED  WITH 
MEASURED  A I R 6 L 0 W CONTOURS 
FROM  FIG.  3-9 . 


during  the  great  geomagnetic  storm  of  May  1967.  The  0V1-10  auroral  zone 
measurements  are  compared  with  solar-wind  velocities  from  Vela  3 and  5 in 
Figures  5-1  and  5-2.  (If  an  auroral  measurement  begins  at  zero,  it  means  that 
an  upper  bound  was  dete- mined.)  There  was  a rather  good  correlation  amono 
the  reported  data,  but  the  data  are  far  from  complete.  Complete  density  data 
for  comparison  with  solar-wind  data  were  provided  by  the  balloon  satellites, 
Explorers  19  and  24,  which  are  also  shown  in  Figures  5-1  and  5-2.  These  data 
were  discussed  i it  Section  5.0. 

The  rates  of  energy  deposition  in  the  auroral  zone  which  Clark  and  Metzger 
deduced  from  airglow  measurements  are  plotted  for  17  orbital  revolutions  of 
the  0V1-10  satellite  during  the  months  February  and  May  1967.  The  ranges  of 
reported  measurements  have  been  extended  by  +20  percent,  to  allow  for  possible 

random  errors  of  measurement.  When  the  solar  wind  is  very  quiet  and  about 

2 

300  km/sec,  the  rate  of  energy  deposition  is  usually  below  5 ergs/cm  sec,  and 
sometimes  below  1 erg/cm  sec.  But  during  disturbed  times,  the  rate  is  higher, 
sometimes  much  higher  (especially  during  the  great  storm  of  May  25  and  26, 
1967).  The  fact  that  the  airglow  as  measured  in  this  experiment  is  much  more 
variable  than  the  solar  wind  illustrates  the  importance  of  detailed  modeling 
of  the  magnetospheric  and  atmospheric  processes. 
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6.0  SUMMARY  AND  CONCLUSIONS 


The  development  of  these  magnetic  field  and  neutral  density  models  has  lead 
us  to  recognize  that  at  present,  because  of  our  lack  of  physical  understanding 
of  many  environmental  features,  any  environmental  model  of  practical  utility 
must  be  semiempirical . The  magnetic  field  model  has  become  a useful  tool  to 
many  investigators  because  it  accurately  represents  many  observed  features  of 
the  magnetospheric  magnetic  field.  We  have  learned  then  that  the  first 
requirement  of  any  environmental  model  is  not  that  it  be  based  on  a complete 
understanding  of  all  the  physical  phenomena  involved,  but  rather  that  it  be 
useful . In  order  to  be  useful  to  groups  whose  hardware  systems  are  influenced 
by  environmental  features  (and  for  most  basic  research  applications),  a model 
of  one  of  those  features  must  possess  several  qualities.  First,  it  must 
quantitatively  and  accurately  represent  the  environmental  feature  it  is  intended 
to  describe.  Second,  it  must  be  analytic.  For  many  applications  it  is  necessary 
that  the  parameter  to  be  described  is  defined  continuously  over  a region  and 
that  its  representation  be  differentiable.  Third,  it  must  be  possible  to 
represent  the  model  with  a computer  code  that  is  fast.  The  model  is  not  useful 
if  its  output  parameters  can  be  obtained  only  after  large  amounts  of  computer 
time  have  been  expended.  This  has  important  implications  for  future  decisions 
on  model  development  and  data  collection.  The  ultimate  goal  of  predicting 
near  earth  environmental  effects  will  not  bo  reached  bv  continuing  to  increase 
the  amount  of  monitored  data  to  the  point  that  effectively  the  whole  environ- 
ment is  monitored  continuously.  Large  data  sets  and  slow  computer  codes  will 
not  prove  useful  in  the  attempt  to  develop  a capability  for  real  time 
monitoring  of  environmental  features  (although  such  codes  and  data  sets  may 
be  of  some  merit  for  increasing  our  understanding  of  some  of  the  physical 
processes  operating  In  the  environment).  Rather,  the  ultimate  predictive 
system  will  operate  with  a small  set  of  well  chosen  input  parameters  and  a set 
of  very  fast  computer  codes.  Fourth , the  model  must  be  usable  over  a large 
region  of  space.  For  example,  models  of  the  upper  atmosphere  and  ionosphere 
should  be  global  while  models  of  the  magnetosphere  should  be  valid  from  the 
subsolar  region  to  the  distant  tail  (at  least  to  lunar  orbit).  Fifth,  the 
model  must  also  be  compatible  with  other  aodels.  This  dictates  that  care  be 
used  in  deciding  on  the  proper  input  and  output  parameters.  It  must  also  be 


62 


constructed  so  that  environmental  data  observed  by  satellite  and  ground 
based  detectors  can  be  used  as  input  parameters. 

The  magnetic  field  modeling  work  has  led  directly  to  two  other  modeling  efforts. 
In  one  the  work  on  magnetospheric  magnetic  fields  (By)  is  being  extended  to 
include  variations  produced  in  By  by  changes  in  the  angle  between  the  geo- 
magnetic dipole  axis  and  the  solar  wind  direction.  (This  angle  influences 
the  gross  shape  of  the  magnetosphere  and  related  magnetic  field  topology.) 

This  work  is  being  done  for  AFWL  with  monitoring  function  provided  by  OSR. 

(AFWL  became  interested  in  the  work  on  magnetic  fields  because  of  their  need 
to  keep  track  of  charged  particles  released  by  nuclear  bursts.  It  was  believed 
that  some  of  their  problems  were  related  to  the  use  of  inaccurate  field  models.) 
This  work  includes  the  modification  of  the  "B,  L"  coordinates  so  that  with  the 
improved  model,  charged  particle  data  can  conveniently  be  organized  out  and 
beyond  synchronous  orbit  (Pfitzer  and  Olson,  1975).  This  work  will  provide  an 
accurate  description  of  the  magnetospheric  magnetic  field  during  "quiet"  magnetic 
conditions  and  provide  a means  for  organizing  charged  particle  data  during  such 
periods.  The  other  effort,  funded  by  ONR,  is  to  construct  a quantitative  model 
of  total  ionospheric  electron  density.  In  both  of  these  efforts  the  model 
philosophy,  described  above,  is  being  used. 

We  believe  that  we  are  now  close  to  having  a first  order  set  of  models  that 
describe  the  near  earth  environment.  The  magnetic  field  and  atmospheric 
density  models  now  exist.  The  procedures  for  developing  the  associated  electric 
field  model  have  been  developed  with  support  from  the  ONR.  These  models 
together  with  the  ionospheric  electron  density  model  now  under  construction 
can  oe  used  to  test  the  possibility  of  predicting  near  earth  environmental 
behavior.  Such  a test  would  be  valuable  in  at  least  two  ways.  First,  it  would 
serve  to  verify  the  accuracy  of  the  models.  Second,  it  would  demonstrate  that 
at  least  some  bulk  parameters  describing  the  environment  can  be  predicted.  The 
second  item  is  particularly  important  since  such  a test  would  call  modeling 
efforts  to  tiie  attention  of  many  people  caring  for  and  optimizing  the  performance 
of  hardware  systems.  One  such  test  of  predictive  capability  might  be  performed 


as  follows.  Data  from  the  SOLRAD-HI  and  HELIOS  satellites  on  solar  wind 
parameters  could  be  used  as  input  to  the  magnetic  field  model.  Variation 
in  magnetic  field  topology  calculated  with  the  model  could  be  used  as  input 
to  the  neutral  atmospheric  density  and  ionospheric  electron  density  models 
to  indicate  the  location  and  extent  of  charged  particle  energy  sources.  The 
satellite  data  would  be  used  directly  to  provide  inputs  on  the  intensity  of 
the  charged  particle  fluxes.  The  ionosphere  and  atmosphere  models  could  then 
be  used  to  "predict"  several  parameters  that  could  be  checked  against  truth 
data  provided  by  polar  orbiting  satellites  such  as  those  used  on  the  DMSP 
satellites.  Such  an  effort  would  at  the  same  time  permit  questions  on  systems 
problems  to  be  answered.  Examples  are:  What  is  the  optimal  data  samoling 
rate?  Is  the  satellite  coverage  adequate?  What  are  the  optimum  orbits  for 
this  monitoring  function?  Are  the  data  the  best  for  this  monitoring/oredicting 
purpose  or  should  new  sensors  be  developed? 

Compatibility  of  models  must  also  be  tested.  Each  model  must  individually  be 
examined  to  guarantee  that  meaningful  input  parameters  have  been  used.  The  set 
of  models  should  also  be  checked  to  make  certain  they  provide  the  proper  output 
parameters.  (The  output  parameters  should  provide  meaningful  input  to  other 
codes  that  will  be  used  to  compute  such  things  as  radio  wave  paths  through  the 
ionosphere,  the  trajectories  of  satellites  and  missiles  moving  in  the  earth's 
atmosphere  and  radiation  dose  rates  for  components  and  man  in  space.) 

In  addition  to  their  practical  applications,  this  zero  order  set  of  models 
should  also  make  possible  several  pure  research  endeavors  that  previously 
could  not  be  considered.  An  example  is  the  study  of  upper  atmospheric  winds 
on  a global  scale.  The  density  of  the  earth's  upper  atmosphere  and  the 
electron  density  in  the  ionosphere  are  necessary  inputs  to  such  a study.  It 
would  be  hoped  that  such  scientific  studies  at  the  basic  research  level  would 
provide  the  new  insights  necessary  to  improve  these  models  and  make  them  even 
more  useful  for  the  prediction  of  near  earth  environmental  effects  on  hard- 
ware systems.  We  are  also  quite  interested  in  the  Interaction  of  electric 
and  magnetic  fields  with  charged  particles  and  would  hope  that  these  models 
could  be  used  to  quantitatively  study  many  charged  particle  phenomena  in  the 
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earth's  magnetosphere.  We  believe  that  the  use  of  quantitative  field  models 
will  make  possible  the  quantitative  explanation  of  several  storm  time  and 
magnetospheric  substorm  associated  phenomena  observed  over  short  periods  of 
time. 

The  realization  that  quantitative  models  are  necessary  for  both  basic  and 
applied  research  led  us  to  propose  that  a conference  be  held  on  "Quantitative 
Magnetospheric  Models".  The  conference  was  held  May  6-8  in  La  Jolla  Shores, 
California,  with  the  American  Geophysical  Union  as  the  prime  scientific 
sponsor.  AFOSR,  ONR,  AIAA,  and  MDAC  also  supported  this  conference.  A 
summary  of  the  conference  is  given  in  Appendix  C.  A summary  of  the  meeting 
was  presented  to  the  10th  ESCAB  symposium  (Vienna),  the  annual  (Washington) 
meeting  of  the  AGU,  and  to  the  IMS  Workshop  at  the  IUGG  (Grenoble)  meeting. 

At  the  IUGG  meeting,  IAGA  formed  a new  working  group  on  Quantitative  Models 
with  W.P.  Olson  as  its  Chairman.  The  models  developed  under  the  present 
contract  will  be  considered  by  that  working  group  for  use  during  the  IMS. 

Since  completion  of  work  on  the  magnetic  field  model,  over  50  requests  for 
information  concerning  the  model  have  been  made  and  over  35  copies  of  the 
computer  subroutine  that  computes  the  magnetic  field  have  been  distributed 
to  users.  In  addition  to  the  successful  tests  of  the  model  listed  above,  the 
model  has  been  used  by  various  groups  for  such  purposes  as  determining  the 
foot  of  the  field  line  that  passes  through  the  ATS  satellite  and  the  location 
of  ground  stations  for  the  GEOS  satellite.  The  model  has  also  been  used  to 
successfully  predict  the  trajectory  of  a rocket  such  that  it  could  emit  a 
beam  of  electrons  and  later  In  Its  flight  be  in  the  proper  position  to  detect 
some  of  those  electrons  moving  back  along  magnetic  field  lines  after  being 
reflected  at  the  mirror  point.  It  is  hoped  that  the  atmospheric  density  model 
will  also  become  a useful  tool  for  the  research  community.  We  firmly  feel 
that  both  models  will  ultimately  have  many  practical  uses. 

Under  this  contract,  a program  to  examine  the  influence  of  the  solar  wind  on 
the  near  earth  environment  was  Initiated.  The  goal  of  ultimately  predicting 
the  behavior  of  the  environment  was  kept  in  mind  throughout  the  performance  of 
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the  contract.  Quantitative  models  of  the  dayside  cusp,  total  maqnetospheric 
magnetic  field,  and  neutral  density  of  the  upper  atmosphere  were  developed. 

We  are  confident  that  ultimately  these  models  or  their  derivatives  will  be 
important  inputs  to  a hardware/software  system  for  quantitatively  predicting 
environmental  behavior  and  effects  on  hardware.  We  have  also  been  nleased  to 
find  that  these  practical  tools  have  already  been  used  by  our  group  and  many 
others  to  perform  interesting  and  vital  basic  research. 

The  investigators  wish  to  thank  AFOSR  for  its  support  of  this  work.  They 
greatful ly  acknowledge  computational  help  provided  by  M.  B.  Baker. 
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7.0  APPENDICES 


7.1  Appendix  A:  SUBROUTINE  BXYZ 

A quantitative  model  of  the  distributed  currents  and  their  associated 
magnetic  field  has  been  developed  (Olson  1974;  Olson  and  Pfitzer,  1974). 

This  model  has  been  used  together  with  a model  of  the  magnetopause  currents 
(Olson,  1969)  and  a dipole  representation  of  the  earth's  main  field  to  re- 
present the  total  magnetospheric  magnetic  field,  For  one  of  a kind 
calculations  (constant  £ contours,  AB  contours,  S^  patterns,  etc.),  £ is 
accurately  determined  by  direct  integration  over  the  current  systems. 

However,  for  repeated  usage  (calculation  of  particle  trajectories,  field 
lines,  drift  shells,  etc.),  direct  integration  is  too  expensive  and  instead 
an  analytic  representation  must  be  used  with  some  loss  of  accuracy. 

The  subroutine  BXYZ  returns  an  analytic  representation  of  $.  The  components 
of  $ are  given  in  series  form.  The  coefficients  to  the  series  were  computed 
by  a generalized  mul ti-variable  orthonormal  fitting  program  (Pfitzer,  1973). 
For  both  the  boundary  and  distributed  currents,  the  field  contributions  were 
obtained  at  over  600  points  in  one  quadrant  of  the  magnetosphere  by  direct 
integration.  (The  model  magnetosphere  considered  here  exhibits  symmetry 
about  the  XZ  and  XY  planes.)  The  fitting  program  was  then  used  to  find  the 
coefficients  of  a series  expansion  for  each  of  the  components  of  $.  The 
average  error  in  the  three  components  of  the  field  (determined  by  comparing 
5 calculated  by  direct  Integration  over  the  two  current  systems  with  the 
values  given  by  the  series  representation)  is  9.8%  or  .63y. 

It  is  necessary  to  use  a large  number  of  terms  In  the  series  expansion  because 
the  field  from  the  distributed  and  magnetopause  currents  is  quite  structured 
especially  in  the  near  earth  region.  The  expansions  are: 
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X,  Y,  and  Z are  the  coordinates  (in  earth  radii)  at  which  the  field  is  to  be 
found.  X is  toward  the  sun,  Z is  perpendicular  to  X and  coincident  with  the 
dipole  (the  model  represents  zero  tilt).  Y is  such  as  to  form  a right-handed 
coordinate  system.  The  constants  D,  E and  F are  adjusted  to  give  the  smallest 
total  error  to  the  least  squares  fit.  R2  = X2  + Y2  + Z2 

The  expansion  is  valid  only  in  the  following  region.  For  points  closer  to 
the  sun  than  X * -5  R^  the  expansion  is  valid  within  a 16  half  sphere 
centered  at  X 3 -5  R^.  For  points  further  down  the  tail  than  -5  R^  the  region 

of  validity  is  defined  by  a cylindrical  surface  about  the  X axis  where  the 

diameter  of  the  cylinder  increases  with  increasing  X.  The  diameter  of  the 
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X must  be  greater  than  -68  R^.  Since  a high  order  series  was  used  the  series 
diverges  rapidly  outside  the  region  of  validity.  The  computer  program  con- 
tains a check  on  the  input  position  where  $ is  to  be  determined.  If  it  is 
outside  the  region  of  validity  the  field  is  set  to  the  dipole  value  and  an 
error  message  is  printed.  The  series  does  not  accurately  represent  the 
external  sources  for  r < 1.5  and  therefore  should  not  be  used  for  deter- 
mining variations  ne^.r  the  earth.  The  external  sources  in  this  region  con- 
stitute less  than  IX  of  the  total  field  and  therefore  the  total  field  is  not 
affected  by  the  loss  of  accuracy  of  the  series. 
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The  subroutine  is  written  in  Fortran  IV.  Its  calling  sequence  is: 

CALL  BXYZ  (X,  B,  BMAG) 

X is  a three-dimensional  input  array  where  X(l),  X(2),  and  X(3)  are 
respectively  the  X,  Y,  and  Z coordinates  (in  R^.)  in  solar  magnetospheric 
coordinates. 

B is  a three-dimensional  output  array  where  B(l),  B(2),  and  B(3)  are 
respectively  the  X,  Y,  and  Z components  of  B.  BMAG  is  the  scalar  magnitude 
of  3.  The  magnetic  field  is  expressed  in  gauss  and  includes  the  contribution 
of  the  boundary  and  distributed  currents  as  well  as  a dipole  main  field. 

Although  the  data  points  input  to  the  least  squares  program  (i.e.,  the  inte- 
gration over  the  currents)  guarantee  that  7 * t =■  0,  the  least  squares 
program  does  not  force  $ • $ to  equal  zero.  Therefore,  $ 0 in  the 
expansion.  It  has  been  verified  thatvx  $ of  the  analytic  series  expansion 
does  indeed  reproduce  approximately  the  currents  which  were  initially  input 
into  the  program. 

If  a better  representation  for  the  main  field  is  required  the  dipole 
representation  can  be  replaced  by  a more  accurate  version  of  the  main  field. 

A listing  of  BXYZ,  a test  program,  and  some  sample  output  follow.  A copy 
of  this  deck  is  available  upon  request. 
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PROGRAM  BTEST!  OUTPUT,  TAPE6  = 0UTPUT  ) 

C THIS  PROGRAM  TESTS  SUBROUTINE  BXYZ 
DIMENSION  X<  3 ),B( 3 > 

LNCNT= 100 
DO  1 1=1,65,6 
X ( 1 )=  1 1 - I 
DO  1 J= 1 , 1 3 , 6 
X(  2 )=  J- 7 
DO  1 K=1 , 13,6 
X(  3 ) = K-7 

C THE  IF-TESTS  DEFINE  THE  REGION  OVER  WHICH  THE  B-F I EL  D EXPANSION  IS 
C DEFINED 

C EXIT  IF  THE  POSITION  IS  OUTSIDE  THE  REGION  OF  VALIDITY 
IF<  X(  1 I.LT.-5 . ) GO  TO  5 

IF< ( ( X< 1 )+5 . )**2  + X( 2 )**2  + X( 3 )**2  I.GT.256  . ) GO  TO  1 
GO  TO  6 

5 IF(((X(2)**2+X(3)**2).GE.((9.*(10.-3.*X(1))/(10.-X(1)))+1.  )** 

*2  . ) . OR . ( X(  1 ) . LT . -68 . ) ) GO  TO  1 
6 CALL  BX YZ( X , B , BM AG  ) 

LNCNT=LNCNT+ 1 

I F ( LNCNT.LT. 55  ) GO  TO  20 

LNCNT=0 

WRITE! 6 , 10  ) 

10  FORMAT!  1 H 1 , / , 5 X , 3H X SM , 7 X , 3H Y SM , 7 X , 3H Z SM , 7 X , 2 HB X , 8 X , 2 HB Y , 8 X , 2 H B Z , 
*8  X , 9 HBM AG , //  ) 

20  WR I TE! 6 , 1 00  ) X , B , BMAG 

100  FORMAT!  1H  ,7FlO  .5  ) 

I CONTINUE 

END 
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SUBROUTINE  BX  YZ( XX, BB, BMAG  ) 

C 

C VERSION  MAY  9,  1973  RECODED  AUGUST  10,  1973 
C CODE  MODIFIED  NOVEMBER  6,  1975 

C THIS  SUBROUTINE  RETURNS  THE  MAGNETIC  FIELD  FROM  THE  MAGNETOPAUSE  CURRENTS, 
C CURRENTS  DISTRIBUTED  THROUGHOUT  THE  MAGNETOSPHERE,  AND  A DIPOLAR 
C REPRESENTATION  OF  THE  MAIN  MAGNETIC  FIELD  OF  THE  EARTH.  A GENERALIZED 
C ORTHONORMAL  LEAST  SQUARES  PROGRAM  DAS  USED  TO  CALCULATE  THE  COEFFICIENTS 
C OF  A POWER  SERIES  THROUGH  SIXTH  ORDER  WITH  EXPONENTIAL  TERMS.  THE 
C POSITION  OF  THE  POINT  WHERE  THE  FIELD  IS  TO  BE  FOUND  IS  INPUT  AS  A 3 
C COMPONENT  ARRAY, XX.  X X ( I ) , X X ( 2 ) , AND  X X ( 3 ) ARE  RESPECTIVELY  THE  X , Y , AND 
C Z COORD  1ANTES  OF  THE  POINT  GIVEN  IN  EARTH  RADII  IN  SOLAR  MAGNETOSPHER I C 
C COORDINATES.  THE  MAGNETIC  FIELD  IS  RETURNED  IN  THE  ARRAY  BB.  THE  X , Y , AND 
C Z COMPONENTS  OF  THE  FIELD  ARE  RESPECTIVELY  BB(  1 ) , BB( 2 ) , AND  BB( 3 1 . THE 
C FIELD  IS  RETURNED  IN  GAUSS.  <X  IS  POSITIVE  IN  THE  SOLAR  DIRECTION.) 

C BM AG  IS  THE  MAGNITUDE  OF  THE  FIELD  IN  GAUSS. 

C 

C THIS  SUBROUTINE  SHOULD  NOT  BE  USED  FOR  X LESS  THAN  -fc8RE.  A TEMPLATE 
C CHECKS  THE  LOCATION  OF  THE  POINT  AND  IF  THE  POINT  IS  OUTSIDE  THE 
C REGION  OF  SERIES  CONVERGENCE  THE  FIELD  IS  SET  TO  THE  DIPOLE  VALUE  AND 
C AN  ERROR  MESSAGE  IS  ISSUED. 

C 

C DO  NOT  USE  THIS  SUBROUTINE  FOR  DETERMINING  VARIATIONS  NEAR  THE 
C EARTH  THE  LEAST  SQUARES  FIT  DOES  NOT  ACCURATELY  REPRESENT  THE 
C EXTERNAL  FIELD  CONTRIBUTION  FOR  R LESS  THAN  1.5  RE.  (THE  EXTERNAL 
C FIELD  IS  LESS  THAN  1 PERCENT  OF  THE  TOTAL  FIELD  IN  THIS  REGION). 

C 

C THE  VALUES  FOR  THE  COEFFICIENTS  IN  THIS  LISTING  REPRESENT  THE  ENTIRE 
C EXTERNAL  FIELD  ( BD  + DS  ) . COEFFICIENTS  REPRESENTING  THE  DISTRIBUTED 
C AND  BOUNDARY  FIELDS  SEPARATELY  ARE  ALSO  AVAILABLE. 

C 

C FOR  REFERENCE  SEE  -THE  MAGNETOSPHER IC  MAGNETIC  FIELD-  BY  W.P.  OLSON 
C AND  K . A . PFITZER  (SEPTEMBER  79  J.G.R),  -THE  DISTRIBUTED  MAGNETOSPHER  1C 
C CURRENT5-  BY  W.P.  OLSON  (SEPTEMBER  79  J.G.R),  AND  -THE  SHAPE  OF  THE 
C TILTED  MAGNETOPAUSE-  BY  W.P.  OLSON  ( J . G .R  .,  79 ,5692,  1 969  ) . 

C 

C PLEASE  ADDRESS  ANY  QU E5T I ON S , C OMMENTS  OR  CRITICISMS  i J P.  OLSON  OR 
C K. A. PFITZER,  MCOONNELL  OOUGLAS  ASTRONAUTICS  COMPANY, 

C 5 30  1 BOLSA  AVE, HUNTINGTON  BE ACH , C AL  I FORN  I A , 9 269  7 
C OR  PHONE  (719)  896-9368  OR  896-3231 
C 

DIMENSION  BB(  3 ),  XX(  3 ),  A(  30  !,B(  30  ),C(  20  ),  D(  20  ),E(  39  ),  F(  39  ) BD  + DS 

DATA  ( At  I ),  1=  1 , 30  )/  7.831  7 9E-0  1 , 2 . 8 80 06 E-0 3 , - 1 . 6 8800 E-05 , BD  + DS 

* 3 . 56869E-08 , -2 . 28690E-03 , 1 . 2969 7E-05  , - 1 . 1 75 7 3E-08  , 3.79779E-06,  BD  + DS 
+-8  .627  00E-09  , 9 . 6 35  8 I E-09  , - 8 . 5 8 65  1 E-0 2 , 2.207  1 3E-09  , 9.9  1626E-07,  BD  + DS 

* 1 . 209 1 9 E-09  , 3.  19819E-07,  2 . 20 8 86 E-0 7 , -6  . 0 7 0 7 0 E-0 3 , 1.92095E-05,  BD  + DS 

+ 6 .96  12  3E-09 , 9.63390E-06  , 2.6  3098E-09  , 1 . 9 7 7 8 3E- 09 , - 1 . 6 8 25  9E-09  , BD  + DS 

* 2.621 1 9 E-0  7 , 1 . 8095  2E-07  ,-2 . 839  19E-06  , 1.995  32E-09,  1.06639E-09,  BD  + D5 

*-2  .97037 E-08,-1 . 38  79  7E-10/  BD+DS 

DATA  (B<  I ),  1=1,30  )/-l  .36002E  + 00,  7 . 8 0269 E-02 , - 1 . 0 08 7 2E -0 3 , BD+DS 

+ 3.881 02E-05 , 9.68926E-02,  9.99081E-09,  3 . 65 6 92 E-05 ,- 3 . 0 889 2 E-0 3 , BD+DS 

* 8.279  18E-05,  2. 26903E-05  , - 2 .01  122E  + 00,  2 . 965  1 3E-02  , - 1 . 5 0385  E-09 , BD  + DS 

*-3  . 76012E-02,  3 . 2 1 39 7 E-09  , - 3 . 7 75  1 6 E-09  , 5 92951E-01,  1.7  7058E-03,  BD  + DS 

*-6  . 37539E-05  , -7 .95695E-03,-l  .02199E-09,  1 . 2 1 1 99 E-09 , -2 . 06  25  2E-0? , BD+DS 
+ 6 . 77865E-05 ,- 1 . 37321 E-03, -2 .95859E-02,  2.93961E-05,  2.10955E-09,  BD+DS 
*-3 . 88689E-09,  2.29031E-09/  BD+DC 

DATA  <C(  I ),  1=  1,20  )/-9 . 36032E-02  , 6 . 9 75 8 1 E -05  , - 2 . 8 9 7 1 9 E-0 8 , BD+DS 
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* 1 . 12935E-09,-2.29679E-07,-6  . 5 1 1 05  E-08  , -9  . 1 185 8 E-0 3,  2. 2025  7E-06, 
*-9 . 7 19  35  E- 1 0 , 5 . 18306E-06,-9 . 8 1855E-09, -2 . 098 12E-09, -3 . 7 8 09  1 E-09, 

* 9.7615  9E-0  9 , 9 . 95 099E-08 , -7 . 85 6 7 1 E-06 , -9 . 8 1 5 1 7 E- 1 1 , 7.93833E-10, 
+ -8 .06  396E-08,-3.28896E-10/ 

DATA  (D<  I ), 1 = 1,20  » / — 1 .69179E+00,  2 . 9 30 38E-02 , -9 . 6 2707E-0H , 

+-5  .02388E-02,  3 . 6 1 602E-09 , -6  .00587E-09,  1 . 1 98 72E-0 1 , - 2 . 7 32  9 3E-0 3 , 

* 3. 1509  1E-05  ,-6 .81062E-03,  2.0091  3E-05  , 8 . 7 1659E-05 , -3 . 932  1 1 E-02 , 

* 6 .29163E-09 ,-l . 79693E-03,-8  . 77802E-03,  2.921  15E-05  , 2 . 09335  E-OH  , 
*-9 .5  725  2E-09 , 1.33622E-09/ 

DATA  (E( I >,  1=1,39  )/  1 . 00 382E+02 , -9 . 5 1 8 30E-0 1 , 2 . 9 7595 E-09  , 

* 1 . 131  97 E-0 7 ,-1  .8  1773E-01,  3 . 9 7 092E-09  , -9 . 1 790 1 E-0  7 , 2 . 6095  7E-09 , 
*-2  . 12733E-07,  1.16209E-07,  1 . 00 7 1 3E  + 0 1 , -9 . 300 35 E-02  , 1.32877E-05, 

* 1 .89977E-09,-l .03666E-02,  1 . 09 399E-05  , -6 . 7 1 7 3 1 E-09 , 1.76089E-05, 
+-1  .92875E-09,  2.5  3223E-09,  9 . 1 26 89E-0 1 , - 1 . 7 1 790E-0 3 , 2.27366E-07, 
+-1  .05062E-09,  1.00609E-07,  3.73228E-07,  8 . 99098E-03 , - 3 . 985 09E-05  , 

* 1 . 23829E-09 , 3.63987E-06,  1.95769E-10,  2.95990E-09,  8.52182E-05, 

•-3.53777E-07,  8.56056E-08,  3 . 3265 7E-07 , - 1 . 92659E-09 , 9.98217E-10, 
*-8 . 1 095  1 E- 1 1 / 

DATA  ( F< I ), 1 = 1 , 39  )/-l . 1 29 3 1 E ♦ 0 2 , - 6 . 95  878E-01  , -2 . 829  8 9E-0  3, 

+ -1  .75919E-05,-3.85656E+00,  2 . 2 1 995 E-02 , - 1 . 19133E-09,  5 .265  31 E-0  2, 
*-5  . 269 19E-09 , -3 . 10689E-09,-l . 0 3980E+0 1 , - 3 . 26 9 1 7 E-02  , 2.55292E-03, 
*-l . 85172E-06,-3.68576E-02,  3 . 7 75  02E-03, -3 . 99  9 33E-05  , 3. 25 3 75 E-09, 
*-6 . 12569E-05 ,-2. 35922E-05 ,-3 .99999E+00,  3. 6 1909 E-02  , 2 .27925 E-09, 

* 9 . 096 10E-0 2,-1 . 6 289 3E-0  3,-9 . 7 37 2 1 E-09 , - 1 . 1 9 1 25 E-02 , -2 . 99  99  8E-0  3, 
*-8 .98886E-08, -9 .99055E-03, -9 . 02203E-05 ,-l . 95  662E-05  , 1 . 77927E-02, 
*-l . 16655E-03,-9 . 795 1 8 E-09 , -2 . 05 600E-0 3 , -3 . 5 98 15 E-05 , -2 . 088 78E-05 , 
*-l  .58959E-09/ 

X = XX<  1 ) 

V=XX( 2 ) 

Z=XX( 3 ) 

V2=V*»2 
Z2=Z**2 
R2=X**2+Y2+Z2 
BX  = 0 . 

BY=0. 

BZ=0. 

C TEST  FOR  LOCATION  OF  INPUT  P05ITI0N 

C IF  THE  LOCATION  IS  OUTSIDE  THE  REGION  OF  SERIES  CONVERGENCE  SET 
C THE  MAGNETIC  FIELD  TO  DIPOLE  AND  PRINT  AN  ERROR  MESSAGE 
C THE  TEMPLATE  WITHIN  WHICH  THE  5ERIES  CONVERGES  IS  GREATER  THAN  THE 
C DEFINED  LOCATION  OF  THE  MAGNETOPAUSE  BOUNDARY 
IF<  X.LT.-5 . ) GO  TO  5 

IF< ( X+5  . )**2*Y2*Z2.GT. 329  . ) GO  TO  50 
GO  TO  6 

5 !F<< Y2  + Z2  ).GE.< 9 .*< 10 .-3. *X  )/(  10  .-X  )+3.  )**2  ) GO  TO  50 

6 EXPXY=EXP(-0.03*R2  ) 

EXPZ  = EXP<-0.015*R2  ) 

11  = 1 
JJ=1 
KK=i 
XB=  I . 

DO  30  1=1,7 
YEXB=XB 
DO  20  J=1 ,9 

irC i + 2»J.ST  10  ) GO  TO  30 


BD  + DS 
BD+DS 
BD+DS 
BD+DS 
BD+DS 
BD  + DS 
BD+DS 
BD  + DS 
BD  + DS 
BD  + DS 
BD  + DS 
BD+DS 
BD+DS 
BD+DS 
BD+DS 
BD+DS 
BD+DS 
BD+DS 
BD  + DS 
BD  + DS 
BD+DS 
BD  + DS 
BD  + DS 
BD+DS 
BD  + DS 
BD+DS 
BD+DS 
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10  BZ- BZ  + < E<  KK  >+F<  KK  )*EXPZ  )*ZEYEXB 

KK=KK+ 1 

I F(  I JK  . G T . 1 1 > GO  TO  20 
BX  = BX  + < A<  I I ) + B ( II  )*EXPXY)*ZEYEXB*Z 
11=11+1 

I F ( IJK  .GT. 10  ) GO  TO  20 

BY=B Y+( C<  JJ  > + D<  JJ  )*EXPXY  )*ZEYEXB*Z*Y 

JJ= JJ+1 

ZEYEXB=ZEYEXB*Z2 
I JK=  I JK  + 2 

I F ( IJK . L E . 12  ) GO  TO  10 
20  YEXB=YEXB*Y2 

30  XB=XB*X 

C ADD  DIPOLE  FIELD  AND  CONVERT  TO  GAUSS 
HO  AR=-31000 . /<  R2*R2*SQRT(  R2  ) ) 

BB(  1 )=<  3 . *X*Z*AR  + BX  )*0 . 00001 
BB( 2 )=<  3-*Y*Z*AR+BY  > + 0.00001 
BB( 3 )=( ( 3 . +Z2-R2  )*AR+BZ  >*0.00001 
BMAG=SORT( BB(l)**2  + BB(2)  + *2  + BB(3>**2> 

RETURN 

C ERROR  EXIT 
50  WRITE! 6,60)  XX 

60  FORMAT  HH  X=  ,E10.3,HH  Y=  ,E10.3,HH  Z=  ,E10.3,76H  IS  OUTSIDE  THE 
♦ VALID  REG  I ON--POWER  SERIES  DIVERGES  BFIELD  IS  SET  TO  DIPOLE) 

GO  TO  HO 
END 
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7.2  Appendix  B:  SUBROUTINE  ATMOS 


A quantitative  global  model  of  the  neutral  density  of  the  upper  atmosphere 
has  been  developed  and  is  discussed  in  the  text  of  this  report.  This  model 
has  been  coded  in  FORTRAN  IV.  It  can  be  used  from  0 km  to  well  above  500  km. 
The  details  and  functions  used  to  construct  the  model  are  discussed  In  the 
body  of  the  text.  The  relationship  between  the  mathematical  symbols  used  in 
the  text  and  the  FORTRAN  names  are  defined  in  Table  B-l. 

The  calling  sequence  for  the  subroutine  is  CALL  ATMOS  (RR,XLAT,XLONG,DATE, 
UT.DEN) 

RR  is  the  altitude  In  km.  The  model  was  constructed  to  primarily 

represent  the  density  between  120-450  km.  Below  120  km  a 7th 
order  fit  to  CIRA  72  is  included  for  completeness.  The  model 
can  be  used  above  450  km;  the  density,  however,  drops  off  more 
rapidly  than  observed. 

XL AT  Is  the  geographic  latitude  In  dearees  (-90.°  to  +90.®) 

+ representing  North. 

XLONG  is  the  geographic  East  longitude  in  degrees  (0.°  to  +360.°) 

DATE  is  the  day  of  year  (1.  to  366.) 

UT  is  the  universal  time  In  hours  (0.  to  24.  hrs) 

3 

DEN  is  the  output  density  in  gm/cm  . 

Additional  parameters  are  available  thru  the  labeled  COMMON  block  PARAM 
COMMON/ P ARAM/ B IGF ,FBAR ,CLZ ,SZ ,E0Z .ALPHA 

These  6 variables  of  the  COtfiON  are  preset  by  a BLOCK  DATA  routine.  Their 
respective  meanings  are  given  In  table  B-l  and  their  default  values  are  given 
in  Table  B-2.  The  default  values  should  be  modified  to  suit  solar  cycle  and 
particle  precipitation  conditions  In  effect  at  the  time  of  the  density 
calculation. 

A listing  of  Subroutine  ATMOS  along  with  a test  program  (ATTST)  and  sample 
output  follow.  A copy  of  this  deck  is  available  upon  request. 
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L 

c 

fO 

a; 

>0 

IO 

c 

3 

o 

o 

03 

(V 

IO 

a> 

3 

■0 

o 

o. 

< 

cr 

LlJ 

a 

r> 

© 

r— 

CO 

2= 

X 

xz 

X 

xz 

<_> 

a 

u 

Table  B-2 

Default  Values  for  Variables  In  COMMON  block  PARAM 


PROGRAM  ATT 51  ( OUTPUT, TAPE6  = OUTPUT  ) 

C THIS  PROGRAM  PROVIDES  A SAMPLE  OUTPUT  FROM  THE  DENS  I T V SUBROUTINE 
COMMON/PARAM/B IGF, F BAR, Cl Z,SZ,EOZ, ALPHA 
DIMENSION  ARRAY!  18  ) 

DATE-O. 

UT=0  . 

ALT=900. 

DO  50  I A = 1 ,2 
ALT=200  + < I A - 1 )*200 
L INE  = 1 3 

10  FORMAT!  1 OH  1 AL  TI  TU  DE=  , F 7 . 2,  7H  DATE= , F5  . 0 , 7H  TIME=,F6.2) 

DO  50  1=1 ,72,3 

IF! L INE . L T . 1 2 ) GO  TO  15 

L I NE=  0 

WRITE! 6,10  ) AL T, D ATE , U T 
16  L iNE=L  INE  + 1 

XLONG  = ( 1-1  )*5 
WRI TE(  6 , 20 ) XLONG 

20  FORMAT!  1HH0E  . L ONG  I TU 0 E= , F 7 . 2 , 5 2H  LATITUDE  VARIES  FROM  +85  TO  -85 
+ IN  STEPS  OF  10  DEG  ) 

DO  30  J=  1 , 1 8 
XL AT  = 95 - J + 1 0 

30  CALL  ATMOS! ALT, XLAT, XLONG, DATE, UT, ARRAY! J ) ) 

WR  I TE ( 6 , 6 0 ) ARRAY 
60  FORMAT! 6E1 1 . 3 ) 

50  CONTINUE 
END 


r>  r>  r>  r>  r> 


SUBROUTINE  ATMOS! RR, XL  AT , XL 0N6 , D ATE , U T , DEN  ) 


i 

* 

i 


C' 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


VERSION  11/25/75 

THIS  DENSITY  MODEL  WAS  DEVELOPED  AS  PART  OF  AIR  FORCE  OFFICE  DF 
SCIENTIFIC  RESEARCH  CONTRACT  NO.  F99620-72-C-0089  . IT  WAS 
CONSTRUCTED  TO  PRIMARILY  DESCRIBE  TOTAL  DENSITY  FOR  THE 
ALTITUDE  RANGE  120-950  KM.  BELOW  120  KM  A 7TH  ORDER  FIT  TO 
CIRA  72  IS  INCLUDED  FOR  COMPLETENESS.  ALTHOUGH  USEABLE  ABDVE 
950  KM,  THE  DENSITY  DROPS  OFF  MORE  RAPIDLY  THAN  OBSERVED. 

THE  DEVELOPEMENT  OF  THIS  FORTRAN  DECK  IS  DESCRIBED  IN  DETAIL 
IN  THE  FINAL  CONTRACT  REPORT  (DEC.  75).  EARLIER  INTERIM 
CONTRACT  REPORTS  DETAIL  THE  DEVELOPEMENT  DF  THE  FUNCTIONS  AND 
DESCRIBE  THE  CORPUSCULAR  ENERGY  SOURCE. 

QUESTIONS  CONCERNING  THE  MODEL  SHOULD  BE  ADDRESSED  TO 
K . A . PFIT2ER  MCDONNEL  DOUGLAS  ASTRONAUTICS  COMPANY  5301  BDLSA 
AVE.  HUNTINGTON  BEACH  CALIFORNIA  92697 
OR  PHONE  (719)  896-3231/896-9368. 

ARGUMENTS 

RR  = ALTITUDE  (IN  KM  ) 

XL  AT  = LATITUDE  (0.  TO  +-90.  DEGREES) 

XLONG  = E.  LONGITUDE  <0.  TO  360.  DEGREES) 

DATE  = DAY  DF  YEAR  (1.  TO  366.  DAYS) 

UT  = TIME  (0.  TO  29.  HRS  OF  UNIVERSAL  TIME) 

DEN  = OUTPUT  DENSITY  <GM/CM**3) 


THE  FOLLOWING  ADDITIONAL  PARAMETERS  ARE  AVAILABLE  THRU 

THE  COMMON  BLOCK  PARAM.  A BLOCK  DATA  ROUTINE  PRESETS 
THESE  PARAMETERS  — THEY  MAY  BE  MODIFIED  TO  SUIT 
EXISTING  ATMOSPHERIC  CONDITIONS 


B IGF 
FBAR 
CLZ 
SZ 

EOZ 

ALPHA 


10  CM  S*t3  L A R FLUX  IN  UNITS  OF  10**-22  WATTS /M**2/HZ 
3 MONTH  AVERAGE  SOLAR  FLUX  SAME  UNITS  AS  BIGF 
IS  THE  NORTHERN  GEOMAGNETIC  CO-LATITUDE  OF  THE  CUSP  CENTER 
IS  THE  HALF-WIDTH  (LONGITUDINAL)  OF  THE  CUSP  HEATING 
REGION 

IS  THE  HALF-WIDTH  (LATITUDINAL)  DF  THE  CUSP  HEATING 
REGION 

IS  THE  PARTICLE  HEATING  PARAMETER 


SEE  APPENDIX  B OF  THE  FINAL  REPORT  FOR  ADDITIONAL  SYMBOL 
DEFINITIONS 


COmmON/PARAM/B IGF, FBAR, CLZ, SZ, EOZ, ALPHA 

DATA  SA1 , SA2,  SA3,  SB  1 , SB  2, SB  3/ . 19  27  ,- . 1 95  7,- . 09  29  , . 009  8 ,-  . 1 3 7 1 , 

*-.0116/ 

D* TA  SRI , SR2, SR3, SR9/- . 02, . 27E-2 , - . 85E-6 , - . 5 9E-9 / 

DATA  CRO .CR1 .CR2 .CR3 .CR9 .CR5 .CR6 .CR7/1 . 2252E-3.- . 8015593E-  1 . 


80 


. i..'i  ,y|ig)iprgi 


IF( ABS( CMLAT-CLX  ).GT.EO ) GO  TO  80 
S=RAD*SZ 

DETERMINE  MAGNETIC  LONGITUDE 

CML0NG=0 . 

IF(C ML AT. LT. 0.00001 )G0  TO  35 
ZC0S  = CT1 *COL ATS*PC 1 C-C OL ATC * ST  1 
ZSIN=CDLATS*SIN( P-Cl ) 

CMLDNG  = ATAN2(  ZSIN,ZCOS ) 

CD  = MAGNETIC  LDNGITUDE  OF  CUSP  CENTER 

CW=P I2*( ( 16  . 6-UT  )/2H  . ) 

IFICW.GT.PI  ) CU=CU-P I 2 

DL AM  = DISTANCE  IN  MAG.  LONG.  FROM  CENTER  OF  CUSP 

DLAM=CMLONG-CU 
I F ( DLAM .LT . 0 . ) DL AM=DL AM+P I 2 
IF( DLAM .GT . P I ) DL AM=P I2-DLAM 
I F ( DLAM .GT. S ) GO  TO  HO 
XL  1 = D1-G8*DLAM/S 
E=CML AT-CLX 
GO  TO  70 
XL  1 = D1-G8 

E=  AC  OS ( CMLATC*COSCLX  + SIN( CMLAT )*SINCL*COS(  DLAM-S ) ) 

CUSP  HEATING  EFFECT 

IF( ABS( E ).LE .EO )C  = ALPHA*XL1*(  1 .+BETA*C0S( PI*E/EO ) )*(  1 .-EXP( ( 120. 
*)/»)) 

TOTAL  DENSITY 

DEN  = 2 . 7E  - 1 1 * V *(  l.+C  )*EXP((  120.— R )/( < ,99+.518*UVF )*SQRT(  R-103.  ) )) 
**<  1 .+SEMIAN  ) 

RETURN 
C DN=  1 

I F ( R . L T . 9 0 . ) CON  = EXP( -(  R-90 . )**2/200.  ) 

CALCULATE  CIRA  ATMOSPHERE 

ALT=((<(((CR7*R+CR6)*R+CR5)*R+CRH>*R+CR3>*R+CR2)*R+CR1)*R 
DEN  = CR0*EXP( ALT  )*(  1 . + SEM I AN*C ON  ) 

RETURN 

END 


•r 


h.  BLOCK  DATA  ATP10SP 

COnnON/PARAn/B IGF,FBAR,CLZ,5Z,EOZ, alpha 

DATA  BlGF,FBAR,CLZ,5Z,EOZ, ALPHA/115. ,115. ,15. ,90. ,2b  . , Ab  / 
END 


* 


j 
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ALT! TUDE=  200. 

00  DATE 

— 0.  TIME 

s 0.00 

E . 

LON6ITUDE= 

0.00 

LATITUDE  vIr 

IES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 85  E-  1 2 

1 78E-.12 

. 1 6^E- 1 2 

.169E-12 

. 179E-12 

. 1 79E-12 

183E-12 

. 187E-12 

. 191E-12 

. 1 95  E- 1 2 

198E-12 

. 202E-12 

. 205  E-  1 2 

207E-12 

. 208  E- 1 2 

. 2 0 8 E - 1 2 

21  IE-12 

. 2 39  E- 1 2 

E . 

LONGITUDES 

15  .00 

LATITUDE  VARIES  FROM  + 85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 186E-12 

. 180E-12 

. 165E-12 

. 167E-12 

. 171E-12 

. 179E-12 

. 1 78E-1 2 

. 181E-12 

. 185E-12 

. 188E-12 

. 192E-12 

. 196E-12 

. 200E- 1 2 

.203E-12 

. 205  E- 1 2 

. 206E- 1 2 

•209E-12 

232E-12 

E . 

LONGITUDES 

30.00 

LATITUDE  VARIES  FROM  -t-85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. I 8 7 E- 1 2 

. 181E-12 

. 165E-12 

. 169E-12 

. 168E-12 

. 171E-12 

. 1 7*+  E-  1 2 

. 1 7 6 E - 1 2 

. 180E-12 

. 183E-12 

. 187E-12 

. 191E-12 

. 196E-12 

.200E-12 

.202E-1 2 

. 2 09  E- 1 2 

. 207E-12 

. 2 3 1 E-  1 2 

E . 

LONGITUDES 

95  .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 188E-12 

. 182E-12 

. 169E-12 

. 16  3E-12 

. 166E-12 

. 168E-12 

. 1 70E-12 

. 173E-12 

. 176E-12 

. 179E-12 

. 183E-12 

. 188E-12 

. 193E-12 

. 1 9 7 E-  1 2 

. 200E-12 

. 209E-12 

. 206E-12 

. 2 3 0 E — i 2 

E . 

LONGITUDES 

60 . 00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 189E-12 

. 181E-12 

. 169E-12 

. 162E-12 

. 165E-12 

. 167E-12 

. 169E-12 

. 1 71E-12 

. 179E-12 

. 177E-12 

. 181E-12 

. 186E-12 

. 191E-12 

. 196E-12 

. 202E- 1 2 

. 2 1 9 E- 1 2 

. 2 1 9 E-  1 2 

.230E-12 

E . 

LDNG ITUDE= 

75  .00 

LATITUDE  VARIES  FROM  +85 

TD  -85  IN 

STEPS  OF 

10 

DEG 

. 191E-12 

. 179E-12 

. 163E-12 

. 163E-12 

. 166E-12 

. 168E-12 

. 170E-12 

. 173E-12 

. 175E-12 

. 179E-12 

. 183E-12 

. 188E-12 

. 192E-12 

. 197E-12 

•217E-12 

.233E-12 

. 229 E- 1 2 

. 2 3 1 E-  1 2 

E. 

LONGITUDES 

90.00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 193E-12 

. 18  IE-12 

. 169E-12 

. 165E-12 

. 169E-12 

. 1 72E-12 

. 175E-12 

. 179E-12 

. 182E-12 

. 185E-12 

•189E-12 

. 199E-12 

. 198E-12 

. 206  E- 1 2 

. 2 35  E — 1 2 

.25  IE-12 

. 232E- 1 2 

.233E-12 

E. 

LONGITUDES 

105 .00 

LATITUDE  VARIES  FROM  +85 

tO  -85  IN 

STEPS  OF 

10 

DEG 

. 195E-12 

. 189E-1 2 

•166E-12 

. 170E-12 

. 175E-12 

. 180E-12 

. 185E-12 

. 189E-12 

. 193E-12 

. 197E-12 

. 20 1 E- 1 2 

209E-12 

207E-12 

. 2 1 7 E- 1 2 

. 25 1 E- 1 2 

. 269E-12 

235E-12 

. 236E-1 2 

t 

LONG  ITUDE  = 

120.00 

LATITUDE  VARIES  FROM  +85 

TP  -85  IN 

STEPS  OF 

10 

DEG 

. 198E-12 

. 189E-12 

. 171E-12 

. 1 77E-12 

. 185E-12 

. 192E-12 

. 198E-12 

. 209  E- 1 2 

. 20 9 E-  1 2 

. 2 1 3 E- 1 2 

.217E-12 

•219E-12 

220E-12 

..230E-12 

.270E-12 

. 286E-  1 2 

. 258E-12 

•239E-12 

E. 

LDNG ITUDE= 

135 . 00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 2 0 ! E — 1 2 

. 196E-12 

. 178E-12 

. 185E-12 

. 195E-12 

205E-12 

. 2 1 9 E- 1 2 

. 2 2 1 E — 1 2 

. 227E- 1 2 

. 2.32 E- 1 2 

. 2 35  E-  1 2 

. 2 36  E- 1 2 

235E-12 

. 2 9 1 E — 1 2 

. 286E-1 2 

. 31  IE-12 

. 289E-12 

. 293E-12 

E . 

LONGITUDES 

150 . 00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

.205E-12 

. 206  E- 1 2 

. 188E-12 

. 199E-12 

. 207E-12 

. 2 1 9 E - 1 2 

230E-12 

. 2 39 E- 1 2 

•296E-12 

.25  IE-12 

. 25  3 E- 1 2 

.253E-12 

.25  IE-12 

250E-12 

•299E-12 

. 332E-12 

. 300E- 1 2 

.298E-12 

E. 

LONGITUDES 

165 .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DFG 

. 208  E-  1 2 

. 2 1 8 E-  1 2 

. 200E- 1 2 

. 202E-1 2 

.21 7E-12 

. 2 32  E-  1 2 

299E-12 

255E-12 

263E-12 

. 268E-12 

. 2 7 1 E-  1 2 

. 2 7 0 E-  l 2 

. 265  E-  1 2 

. 259E-12 

293E-12 

. 331E-12 

. 308E-12 

. 259E-1  2 

84 


ALTITUDE2  200.00  DA 
E.  LONGITUDE2  180.00 


.21 1 E- 1 2 

.231 E— 1 2 

.256E-12 

268E-12 

. 2 7 7 E- 1 2 

269E-12 

E . 

LONGITUDE2 

195 . 00 

212E-12 

299E-12 

. 265  E- 1 2 

. 2 7 7 E — 1 2 

. 286  E- 1 2 

276E-12 

E . 

LONGITUDE2 

210.00 

212E-12 

. 255E-12 

. 268E- 1 2 

.281E-12 

. 290E-1  2 

. 279E-12 

E . 

LONG  1 TUOE= 

225 .00 

.210E-12 

.261E-12 

. 268  E- 1 2 

.28  IE-12 

.289E-12 

.279E-12 

E . 

LONGITUDE2 

290.00 

. 207E- 1 2 

. 25  OE- 1 2 

263E-12 

. 2 75  E- 1 2 

. 289E- 1 2 

. 2 75  E- 1 2 

E . 

LONGITUDE2 

255  .00 

201E-12 

239E-12 

.255E-12 

.266E-12 

.276E-12 

. 26  7 E- 1 2 

E. 

LONGITUDE2 

270.00 

. 196E-12 

•213E-12 

. 293E-12 

. 25  9 E- 1 2 

. 269E- 1 2 

. 25  7 E- 1 2 

E. 

LONGITUDE2 

285  .00 

. 1 90E-I2 

. 192E-12 

231E-12 

290E-12 

.252E-12 

297E-12 

E. 

LONG 1TU0E= 

300.00 

. 185E-12 

. 183E-12 

.218E-12 

. 226  E- 1 2 

.290E-12 

•236E-12 

E . 

LONGITUDE2 

315.00 

. 183E-12 

. 173E-12 

.207E-12 

. 219E-12 

.228E-12 

. 2 2 6 E — 1 2 

E. 

LONGITUDE2 

330.00 

. 183E-12 

. 1 79E-12 

. 197E-12 

•203E-12 

. 219E-1 2 

-218E-12 

E . 

LONG ITUDE= 

395  .00 

. 189E-12 

. 176E-12 

. 189E-12 

. 199E-12 

. 211E-I2 

. 2 1 2E- 1 2 

= 0.  TIME2  0.00 

LATITUDE  VARIES  FROM  *65 
. 2 15  E- 1 2 . 209E-  1 2 

.277E-12  .283E-12 

. 280E- 1 2 . 3 1 7 E-  1 2 

LATITUDE  VARIES  FROM  +85 
. 239E- 1 2 .216E-12 

.287E-12  . 292E-  1 2 

. 2 7 1 E- 1 2 . 30 1 E-  1 2 

LATITUDE  VARIES  FRDM  +85 
.252E-12  . 226E- 1 2 

.291E-12  . 297E-  1 2 

.267E-12  .289E-12 

LATITUDE  VARIES  FROM  *85 
. 26  1 E- 1 2 . 2 36 E-  1 2 

.290E-12  . 296E-  1 2 

. 265  E- 1 2 . 269E- 1 2 

LATITUDE  VARIES  FRDM  *85 
.262E-12  .299E-12 

. 285  E- 1 2 .290E-12 

. 262E- 1 2 . 25  7 E-  1 2 

LATITUDE  VARIES  FROM  +85 
. 259E-  I 2 . 2*4 5 E- 1 2 

•275E-12  .280E-12 

.256E-12  298E-12 

LATITUDE  VARIES  FROM  *85 
•238E-12  .237E-12 

■262E-12  .267E-12 

. 298 E- 1 2 .239E-12 

LATITUDE  VARIES  FROM  *85 
.223E-12  . 226E- 1 2 ( 

298E-1  2 . 25  3E- 12 

239E-12  . 2 3 1 E-  I 2 

LATITUDE  VARIES  FROM  *85 
. 206  E- 1 2 . 2 06  E-  1 2 

• 233E-12  . 238E-  1 2 

. 23 1 E- 1 2 .229E-12 

LATITUDE  VARIES  FROM  *85 
. 1 85  E- 1 2 . 1 8 5 E-  1 2 

• 219E-12  . 229E-  1 2 

•223E-12  .2I8E-12 

LATITUDE  VARIES  FRDM  +85 
. I TOE- 1 2 -176E-12 

. 208E- \ 2 .212E-12 

-217E-I2  . 2 1 9 E-  1 2 

LATITUDE  VARIES  FROM  +85 
. 166E- 12  .172E-12 

-198E-I2  . 202E- 1 2 

. 2 1 2 E — 1 2 . 2 1 OE - 1 2 


TO  -85  IN 

STEPS  OF 

10 

DEG 

226E-12 

. 292E-1 2 

285E-12 

.283E-12 

208E-12 

.258E-12 

TO  -85  IN 

STEPS  DF 

10 

DEG 

232E-12 

. 299E- 1 2 

299E-12 

. 292E-1  2 

309E-12 

. 26  OE-  1 2 

TO  -85  IN 

STEPS  OF 

10 

DEG 

235E-12 

253E-12 

299E-12 

.296E-1 2 

295E-12 

. 26 1 E- 1 2 

TD  -85  IN 

STEPS  OF 

10 

DEG 

235E-12 

.252E-12 

298E-12 

.296E-12 

283E-12 

. 260E- 1 2 

TO  -85  IN 

STEPS  OF 

10 

DEG 

232E-12 

. 298E- 1 2 

292E-12 

. 290E- 1 2 

272E-12 

.259E-12 

TD  -85  IN 

STEPS  OF 

10 

DEG 

229E-12 

. 29 1 E-  1 2 

282E-12 

. 28 1 E-  1 2 

261E-12 

.257E-12 

TO  -85  IN 

STEPS  OF 

10 

DEG 

223E-12 

. 2 3 1 E-  1 2 

270E-12 

.269E-I 2 

253E-12 

.259E-12 

TO  -85  IN 

STEPS  OF 

10 

DEG 

219E-12 

220E-12 

255E-12 

. 25  5 E-  1 2 

293E-12 

. 25  OE-  1 2 

TO  -85  IN 

STEPS  OF 

10 

DEG 

200E-12 

. 209E-1  2 

290E-12 

. 29 1 E-  1 2 

235E-12 

.297E-12 

TD  -85  IN 

STEPS  OF 

10 

DEG 

191E-12 

. 199E-12 

227E-12 

. 228E-1 2 

227E-12 

.293E-12 

TD  -85  IN 

STEPS  DF 

10 

DEG 

189E-12 

. 191E-12 

. 2 1 5 E — 1 2 

.21 7E-12 

220E-12 

. 239E-1 2 

TO  -85  IN 

STEPS  OF 

10 

DEG 

1 78E-12 

. 189E-12 

206E-12 

. 209E- 1 2 

2 1 5 E- 1 2 

. 2 36  E-  1 2 

85 


*LTITUDE=  *400.00  0ATE=  0.  TIME=  0.00 


E. 

LONGITUDE: 

0.00 

LATITUDE  VARIES  EROH  ->85 

TO  -85  IN 

STEPS  OF 

10 

OEG 

.107E-IM 

. 1 06E- 1*4 

.993E-I5 

. I 05 E-I *4 

1 1 IE-1‘4 

. 1 1 7 E-  1 *4 

- 1 23E-1H 

. 1 29E-  1*4 

. I 3*4E-I*4 

. I *4 OE-  1 *4 

1 *45  E- I *4 

. 1 5 OE- 1 *4 

. 155E-  1*4 

. 1 5 8 E-  1 *4 

-160E— 1*4 

. 1 5 9 E — 1 *4 

1 62E- 1 *4 

. 18  IE-1‘4 

E. 

LONGITUDES 

15  .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. I 07E-  1*4 

. 1 06E-  1*4 

•983E-I5 

. I 0 1 E- 1 *4 

. 107  E—  I *4 

. 1 I 2E- 1 *4 

. 1 1 6 E- 1*4 

. 1 2 1 E- 1 *4 

. 125E-1*4 

. 1 3 1 E- 1*4 

. I 36  E- I *4 

. 1 *4  2E- 1 *4 

. 1 *4  7E- 1 *4 

. 1 5 2 E- 1 *4 

- 155E— 1*4 

. 1 5 6E- 1 *4 

1 5 9E- 1 *4 

. 1 79E-  1*4 

E . 

LONGITUDES 

30.00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 08  E-  1 *4 

. 1 07  E- 1 *4 

. 9 7*4  E- 1 5 

. 989E-I 5 

, 1 03E-I  *4  • 

. 1 0 7 E- 1 *4 

. 1 1 OE- 1 *4 

. 1 1 *4 E- 1 *4 

. 1 1 8 E-I  *4 

. 1 2 3 E - 1 *4 

. 1 29E-  1 *4 

. 1 35  E-  1 *4 

. 1 *4 1 E- 1 *4 

. I *4  7E-  1 *4 

. 1 5 1 E-  1*4 

. 1 5 3E- 1 *4 

. 156E-1*4 

. 1 7 7 E-  1 *4 

E. 

LONGITUDES 

*45  .00 

LATITUDE  VARIES  FROM  *85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 0 8 E- 1 *4 

. 107E-  1*4 

. 965  E-  1 5 

•970E-15 

. 1 OOE- 1 *4 

. 1 0 3 E- 1 *4 

. 1 06E-  1*4 

. 1 09E-  1*4 

. 1 1 3E-  1*4 

. 1 1 8E- 1 *4 

. 1 2 3 E — 1 *4 

. 1 30E-  1*4 

. 1 3 7 E- 1 *4 

. 1*4  3E-  1 *4 

. 1 *4  8 E — 1*4 

. 1 5 2 E- 1 *4 

. I 5 5 E- 1*4 

. 1 7 7 E-  1 *4 

E . 

LONGITUDES 

60.00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 09E-  1 *4 

. 1 06  E-  1 *4 

•957E-15 

. 96  1 E- 1 5 

•993E-15 

. 1 02  E- 1 *4 

. 1 0*4  E-  1 *4 

. 1 07E- 1*4 

. 1 1 IE-1‘4 

. 1 1 5 E- 1 *4 

. 121  E-l  *4 

. 1 28E- 1*4 

. 1 35  E-  1 *4 

. 1*4  IE-1‘4 

. 1 *4  9E-  1 *4 

. H2E-1*4 

. 1 6 2E- 1 *4 

. 1 7 7 E- 1 *4 

E . 

LONGITUDES 

75  .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 1 OE-  1*4 

. 1 05  E-  1*4 

. 95  *4 E—  1 5 

. 969E-  1 5 

. I OOE-l  *4 

. 103E-1M 

. 106E-1*4 

. 109E-1*4 

. 1 1 3 E-  1 *4 

. 1 I7E-1*4 

. 123E-19 

. 1 30E-  1 *4 

. 1 36  E- 1 *4 

. 1*4  3E-  1*4 

. 1 6*4E-  1 *4 

. 1 80E-1  *4 

. 1 7 1 E- 1 *4 

. 178E-19 

E. 

LONGITUDES 

90.00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 1 2E- 1 *4 

. 1 06E-  1 *4 

•967E-15 

. I OOE-I  *4 

. 1 05  E-  1 *4 

. 109E-19 

. 1 1 3 E- 1 *4 

. 1 1 7 E-  1 *4 

. 1 2 1 E- I *4 

. 1 26E- 1 *4 

. 1 32E- 1 *4 

. 1 3 8 E- 1 *4 

. 1 *4*4E-  I *4 

. 1 5*4  E- I *4 

. 1 8 ?E-  1 *4 

. 198E-1*4 

. 1 80E-  1*4 

. 1 7 9E-  1 *4 

E . 

LONGITUDES 

105 .00 

LATITUDE  VARIES  FROM  -^85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 1 *4E-  1*4 

. 1 1 OE- 1 *4 

. 1 OOE-  1 *4 

. 1 06E- 1 *4 

. 1 1 3 E-  1 *4 

. 1 20E-1  *4 

. 1 26E-  1 *4 

. 132E-1*4 

. 1 38  E- 1 *4 

. 1 *4  3E-I  *4 

. 1 *4  9E- 1 *4 

. 1 5 *4  E-  1 *4 

. 1 5 8E-  1*4 

. 1 6 9E- 1 *4 

. 202E-  1 *4 

. 2 1 *4  E — l *4 

. 18*4E-1*4 

. 1 8 2 E — 1 *4 

E. 

LONGITUDES 

120.00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 1 6 E — 1 *4 

. 1 15E-1*4 

. 1 06E-  1*4 

. 1 1 5 E- 1 *4 

. 1 25  E- 1 *4 

. 1 35  E- 1 *4 

. 1 *4  5 E-  1 *4 

. 1 5*4 E-  1 H 

. 1 6 2E-  1*4 

. I 6 8 E-  1 *4 

. 1 7*4E-1*4 

. 1 7 7 E-  1 *4 

. 1 7 9E-  1 *4 

. 189E-1H 

. 228E-  1*4 

. 2*4 1 E- 1 <4 

. 2 I 0 E — 1 *4 

. 1 86  E- 1 *4 

E . 

LONGITUDES 

135  .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 19E-1*4 

. 1 22E-  1*4 

. 1 1 5 E — 1 *4 

. I 26 E-  1 *4 

. 1*4  IE-1‘4 

. 155E-19 

. 1 69E-  1 *4 

. 182E-  1*4 

. 1 92E-  I *4 

. 200E-  1 *4 

. 2 05  E-  1 *4 

. 20 7 E- 1 *4 

. 205 E-  1 *4 

. 2 1 OE-  1*4 

. 25  3 E — I *4 

. 2 7 3E  — 1 *4 

. 2*4  3E-  1 *4 

. 1 9 1 E-  1 *4 

E . 

LONGITUDES 

150 .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 2 2 E-  1 *4 

. 1 32E-  1*4 

. I26E-I*4 

. I 38  E- 1 *4 

. 158E-1H 

. 1 7 7 E-  1 *4 

. 1 96  E- 1 *4 

. 2 1 2E-  1*4 

. 226E- 1 *♦ 

. 2 35  E-I  *4 

. 2*4  OE-  1 *4 

. 239E-  1*4 

. 2 3*4 E-  1 *4 

. 230E-  1*4 

. 27  3E- 1 *4 

. 3 0*4  E— I *4 

. 25  9E-  1 *4 

. 1 9 7 E-  1 *4 

E . 

LONGITUDES 

165 .00 

LATITUDE  VARIES  FROM  +85 

TO  -85  IN 

STEPS  OF 

10 

DEG 

. 1 25  E- 1*4 

, 1 *4*4  E-  I *4 

. 139E-!*4 

. 15  1E-I*4 

. I 75  E- 1*4 

. 1 99E-  1 *4 

. 222E-  1*4 

. 2*4  3 E-  1 *4 

-259E-1*4 

. 269E- 1 *4 

. 2 7*4  E- 1 *4 

. 27  2E- 1*4 

. 25  3E- 1 *4 

. 2H9E-  1 *4 

. 2 8 1 E- I *4 

. 3 1 2E- 1 *4 

. 2 7 I E- 1 *4 

. 20*4 E- 1 *, 
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*LTITUDE=  MOO . 00 


E. 

LONG ITUDE= 

180.00 

. 128E-1M 

. 156E-1N 

.2M5E-1M 

.269E-1M 

.288E-1M 

. 270E-1M 

E. 

LONG ITUOEs 

195.00 

. 129E-1M 

.168E-1N 

.262E-1M 

.288E-1M 

•306E-1M 

.285E-1M 

E. 

LONGITUDES 

210.00 

. 130E-1M 

. 177E-1M 

.270E-1M 

.297E-1M 

•315E-1M 

.292E-1M 

E. 

LONGITUDES 

225 .00 

. 128E-1M 

. 182E-1M 

. 268E-1N 

.295E-1M 

. 313E-1M 

.291E-1M 

E. 

LONGITUDES 

2M0. 00 

. 125E-1M 

.173E-1M 

.258E-1M 

•28ME-1M 

. 303E-1M 

.282E-1M 

E. 

LONGITUDES 

255.00 

. 121E-1M 

. 158E-1M 

.2M2E-1M 

.265E-1M 

.28ME-1M 

■267E-1M 

E. 

LONG ITUOEs 

270.00 

.117E-1N 

. 1M0E-1M 

.221E-1M 

.2M1E-1M 

.261E-1M 

■2M7E-1M 

E. 

LONG ITUOEs 

285.00 

. U2E-1M 

. 121E-1M 

. 198E-1N 

•215E-1M 

■237E-1M 

■227E-1M 

E. 

LONGITUDES 

300.00 

. 108E-1M 

. 1 13E-1M 

. I 77E-1M 

. 190E-1M 

.21 3E-1M 

.207E-1M 

E . 

LONG ITUOEs 

315 .00 

. 106E-1H 

. 10ME-1M 

. 158E-1M 

. 169E-1M 

. 193E-1M 

. 1 90E- 1 m 

E. 

LONGITUDES 

330.00 

. 105E-1M 

. 10ME-1M 

. 1M3E-1M 

. 152E-1M 

. 1 77E-1M 

. 177E-1M 

E . 

LONG ITUDE= 

3M5 .00 

. 106E-1M 

. 105E-1M 

. 1 32E-1M 

. 139E-1M 

. 16ME-1N 

. 166E-1M 

= 0.  TIME=  0.00 

LATITUDE  VARIES  FROM  +65 
. 155E-1N  .161E-1H 

.28TE-1M  .299E-1N 

.2TME-1M  .30ME-1M 

LATITUOE  VARIES  FR0M--+85 
. 1 75E-1M  .1T1E-1H 

. 308E-1M  .321E-1M 

.268E-1N  .289E-1M 

LATITUOE  VARIES  FROM  +85 
.192E-1M  .182E-1M 

.318E-1M  .331E-1M 

.265E-1M  .272E-1M 

LATITUDE  VARIES  FROM  *85 
.201E-1M  .193E-1H 

•316E-1N  .329E-1M 

.263E-1N  .25ME-1M 

LATITUDE  VARIES  FROM  +85 
.199E-1M  .198E-1M 

•30ME-1M  .3I6E-1M 

•256E-1N  .239E-1M 

LATITUDE  VARIES  FROM  *85 
.189E-1M  .195E-1M 

.283E-1M  . 295E-1H 

.2MME-1M  .22ME-IM 

LATITUDE  VARIES  FROM  +85 
.171E-1M  .183E-1M 

•257E-1M  .267E-1H 

. 229E-1M  .211E-1M 

> 

LATITUDE  VARIES  FROM  +85 
. 155E-1M  .168E-1M 

. 228E-1M  . 238E-1M 

•213E-1M  .197E-1M 

LATITUDE  VARIES  FROM  +85 
. 1 38E-1M  .IM5E-IM 
.202E-1M  .210E-1M 

. 198E-1M  ..186E-1M 

LATITUDE  VARIES  FROM  +85 
•119E-1N  . 12ME-1M 

.178E-1M  . 186E-1H 

. 185E-1M  .176E-1M 

LATITUDE  VARIES  FROM  +85 
. 105E-1M  . 1 1ME-1M 

. 160E-1M  .166E-1M 

. 1 7ME-1M  .169E-1H 

LATITUDE  VARIES  FROM  +85 
•10IE-1M  . 1 09E-IM 

. 1M5E-1M  .151E-1M 

•I66E-1M  .16ME-1M 


TO  -85  IN 

STEPS  OF 

10 

DEG 

190E-1M 

.218E-1M 

30ME-1M 

.300E-1M 

276E-1M 

.209E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

200E-1M 

. 232E-1M 

325E-1M 

. 320E-1M 

27HE-1H 

•212E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

205E-1M 

•238E-1M 

335E-1M 

. 329E-1M 

266E-1M 

. 21 2E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

20ME-1M 

.237E-1M 

333E-1M 

. 328E-1M 

25ME-1M 

•212E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

199E-1M 

.229E-1M 

321E-1M 

.316E-1M 

2M0E-1M 

.21  1E-1M 

TD  -85  IN 

STEPS  OF 

10 

DEG 

191E-1M 

.215E-1M 

299E-1M 

. 295E-1M 

226E-1M 

.208E-1M 

TD  -85  IN 

STEPS  OF 

10 

DEG 

180E-1M 

. 198E-1M 

272E-1M 

.270E-1M 

21ME-1M 

.20ME-1M 

TO  -85  IN 

STEP 5 OF 

10 

DEG 

165E-1M 

. 1 79E-1M 

2M3E-.M 

.2M2E-1M 

202E-1M 

. 199E-1M 

TD  -85  IN 

STEPS  OF 

10 

DEG 

1M7E-1M 

. 162E-1M 

215E-1M 

.2 16E-1M 

191E-1M 

. 195E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

13ME-1M 

. 1M6E-1M 

191E-1M 

. 193E-1M 

180E-1M 

. 190E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

12ME-1M 

. 1 3ME-1M 

1 71E-1M 

. 175E-1M 

172E-1M 

. 186E-1M 

TO  -85  IN 

STEPS  OF 

10 

DEG 

1 1 7E-1M 

. 12ME-1M 

156E-1M 

. 16 1E-1M 

l ' 6 E-  1 M 

. 183E-1M 
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7.3  Appendix  C:  Summary  of  the  La  Jolla  Conference  on  Quantitative 
Magnetospheric  Models 

Introduction 


The  meeting  was  held  in  La  Jolla  Shores,  California,  from  6-8  May  1975.  There 
were  five  sponsors:  the  American  Geophysical  Union,  The  American  Institute  of 
Aeronautics  and  Astronautics,  Air  Force  Office  of  Scientific  Research,  the 
Office  of  Naval  Research,  and  McDonnell  Douglas  Astronautics  Company.  The 
meeting  was  divided  roughly  into  three  parts,  with  one  day  devoted  to  each. 

Further,  each  day  was  divided  such  that  approximately  two- thirds  of  the  time 
was  devoted  to  presentation  of  papers,  while  one-third  of  the  time  was  used  for 
discussions  and  working  group  activities.  The  first  day  was  devoted  to  magnetic 
fields,  the  second  day  to  electric  fields,  and  the  last  day  to  charged  particles. 
This  summary  will  roughly  follow  the  order  of  the  meeting.  Many  of  the  comments 
made  concerning  the  magnetic  field  apply  equally  well  to  the  electric  field. 

There  were  85  registered  attendees  v/ith  six  foreign  participants.  One  of  the  pur- 
poses of  the  meeting  was  to  bring  together  those  people  working  on  the  actual 
development  of  models  with  a cross-section  of  those  workers  collecting^observa- 
tional  data.  The  data  is  necessary  for  the  proper  development  of  empirical  models. 
There  were  also  several  model  "users"  present.  Their  current  and  future  needs 
make  it  possible  to  test  existing  models  and  to  suggest  which  future  models  will 
be  most  valuable.  The  meeting  discussions  were  limited  to  quanti tati ve  models 
and  associated  observational  data  sets.  A quantitative  model  was  defined  as 
one  that  yields  numerical  information  which  agrees  well  with  quantitative  obser- 
vations of  a magnetospheric  feature. 

Magnetic  Fields 

In  contrast  to  the  electric  field,  the  magnetospheric  magnetic  field  is  relatively 
well-behaved  and  its  variability  fairly  well  understood.  Modelling  of  magneto- 
spheric  magnetic  fields  is  well  ahead  of  modelling  of  the  magnetospheric  electric 
field  and  various  particle  structures  in  the  magnetosphere.  This  is,  in  part, 
because  the  magnetic  field  is  measured  directly  (electric  fields  are  still 
typically  inferred  from  particle  motions).  Also  the  sources  of  the  magneto- 
spheric magnetic  field,  the  several  magnetospheric  current  systems,  are  fairly 
well  known  and  understood.  The  most  Important  are:  the  magnetopause  currents, 
formed  by  the  deflection  of  shocked  solar  wind  particles  off  of  the  magnetospheric 
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magnetic  field,  the  tail  current  system, produced  mainly  by  charged  particles 
drifting  across  the  plasma  sheet  region  of  the  tail,  and  the  ring  current 
system  which  is  formed  by  trapped  particles  drifting  around  the  earth. 

Sugiura  has  suggested  that  at  least  most  of  the  time,  the  quiet  time  ring 
current  system  is  contiguous  with  the  tail  current  system.  There  are  many  data 
sets  that  have  been  used  to  develop  the  current  models  of  the  magnetospheric 
magnetic  field,  3.  Sugiura,  Hedgecock,  McPherron,  and  Cain  talked  about  obser- 
vations of  it.  Cain  discussed  the  modelling  of  the  earth's  main  field  and  its 
importance  for  magnetospheric  phenomena  as  was  described  at  the  Zmuda  Memorial 
Conference.  A key  point  for  magnetospheric  physicists  was  that  the  present 
International  Geomagnetic  Reference  Field  currently  only  models  "structures" 
with  longitudinal  extent  greater  than  45°.  This  is  probably  inadequate  for  the 
proper  quantitative  study  of  high  latitude  magnetospheric  (and  ionospheric) 
phenomena.  Sugiura  discussed  the  procedure  of  subtracting  the  magnitude  of  the 
earth's  main  field  from  the  magnitude  of  the  observed  field  to  obtain  AB.  It  has 
become  an  important  tool  in  understanding  the  magnetic  effects  of  currents  flowing 
in  the  inner  magnetosphere.  Magnetic  observations  at  synchronous  orbit  also 
provide  excellent  tests  of  the  models  in  that  region  of  the  magnetosphere. 
McPherron  pointed  out  that  the  variations  seen  at  synchronous  orbit  are  quite 
sensitive  to  geomagnetic  longitude.  The  analysis  of  the  HEOS  data  by  Hedgecock 
and  his  colleagues  has  shown  that  in  the  high  latitude,  high  altitude  magneto- 
sphere,, the  magnetopause  is  more  flared  than  predicted  by  the  models  whose  shape 
was  determined  using  the  pressure  balance  condition  (specular  reflection).  This 
data  also  suggests  that  although  charged  particles  may  enter  the  magnetosphere  at 
high  latitudes  over  a cleft  or  cusp-like  region,  the  magnetic  field  topology  in 
that  region  is  nr  like  a funnel  in  shape  and  does  not  extend  more  than  two 
hours  on  either  side  of  the  noon  magnetic  meridian.  Hedgecock's  group  has  also 
initiated  the  use  of  a "hinged"  coordinate  system  that  near  earth  is  like  an 
ordinary  magnetic  coordinate  system,  while  in  different  regions  of  the  magneto- 
sphere it  becomes  either  a solar  magnetic,  or  solar  magnetospheric  coordinate 
system.  With  the  use  of  this  coordinate  system,  it  was  possible  to  do  a much 
better  job  of  organizing  the  magnetic  field  data.  (This  coordinate  system  takes 
into  account  "tilt"  effects.  Tilt  refers  to  the  varying  angle  between  the  geo- 
magnetic dipole  axis  and  the  solar  wind  direction.) 
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Several  models  of  ft  were  presented  by  their  authors  and  reviewed  by  Wa-ker. 

These  included  the  models  of  Sugiura  and  Poros,  Mead  and  Fairfield,  Choe  and 
Beard  and  Sullivan,  Olson  and  Pfitzer,  Willis  and  Pratt,  Alekeev  and  Shabansky, 

and  Voigt.  Most  of  these  models  are  discussed  in  some  detail  in  the  accom- 

panying extended  abstracts.  A general  set  of  requirements  that  an  all  purpose 
model  should  have,  include;  the  magnetopause,  tail  and  quiet  time  ring  currents 
must  be  included  as  must  the  effects  of  tilt.  The  model  must  be  good  through- 

out  a large  portion  of  the  magnetosphere.  It  must  be  capable  of  predicting  ii 

during  quiet  to  moderately  disturbed  conditions.  None  of  the  models  of  & devel- 
oped to  date  meet  all  of  these  requirements. 

In  addition  to  these  general  requirements,  there  are  some  more  specific  details 
that  should  probably  be  included  in  any  future  modelling  efforts.  Potemra 
pointed  out  the  importance  of  field  aligned  currents  in  the  high  latitude  region. 
Although  they  are  presently  not  well  understood  and  their  coupling  of  the  magneto- 
spheric  current  systems  not  presently  known,  it  is  clear  that  at  high  latitudes 
they  exert  some  influence  on  the  structure  of  the  magnetic  field.  Likewise,  at 
high  latitudes,  Hr.dgecock  shaved  that  the  observed  magnetospheri c boundary  is 
more  distant  than  that  given  by  the  older  models.  Either  a better  theory  for 
the  shape  of  the  magnetopause  should  be  developed  or  these  observations  be  input 
to  an  empirical  model.  It  is  also  becoming  increasingly  obvious  that  several 
magnetospheri c phenomena  are  directly  linked  to  the  interplanetary  magnetic  field. 
It  is  therefore  appropriate  to  develop  a model  that  at  once  includes  both  the 
magnetospheri c and  interplanetary  fields.  However,  the  difficulties  with  the 
mathematical  representation  of  such  a field  are  formidable.  Another  requirement 
of  these  models  is  that  they  be  differentiable.  This  is  necessary  if  they  are 
to  be  used  to  interpret  the  drift  and  energization  of  charged  particles.  It  was 
suggested  that  specialized  models  might  be  developed  for  various  local  regions 
of  the  magnetosphere.  However,  such  models  cannot  be  used  conveniently  with 
other  models.  Since  adding  them  together  and  truncating  outside  their  region 
of  validity  causes  severe  problems  with  their  differentiability.  The  same  is 
true  of  trying  to  represent  the  interplanetary  field  and  B simultaneously.  It 
was  agreed  that  any  future  model  should  be  "global",  differential'10,  and  include 
at  least  the  magnetopause,  tail,  and  quiet  time  ring  current  system  contribution, 
and  the  ti It  effect. 


90 


There  were  some  talks  given  on  model  uses.  Seely  reported  that  the  Stanford 
group  has  used  the  Olson-Pfitzer  model  to  locate  the  crossing  of  whistler  ducts 
in  the  equatorial  plane.  McPherron  si.uwed  the  importance  of  including  both  the 
magnetic  effects  of  the  trapped  currents  and  the  effects  of  the  tilt  of  the  geo- 
magnetic axis  to  the  solar  wind  direction  for  the  proper  description  of  the  vari- 
ations  In  13  seen  by  ATS- 1 and  ATS-6  at  synchronous  orbit.  The  Mead-Fai rf ield 
model  accurately  represents  dai ly  variations  in  the  observed  field  perpendicular 
to  the  northward  direction  which  is  most  sensitive  to  tilt.  Hendrickson  told 
of  several  model  tests  provided  by  an  electron  echo  experiment  in  which  electrons 
were  accelerated  away  from  a rocket  and  later,  mirroring  along  magnetic  field 
lines,  detected  by  the  same  rocket.  In  order  to  observe  the  "echo"  it  was 
necessary  to  accurately  predict  both  the  bounce  time  and  the  drift  rate  of  the 
electrons.  The  Olson-Pfitzer  mode!  worked  well  for  this  purpose. 

A general  conclusion  on  the  magnetic  field  modelling  was  that  further  advances 
should  be  directed  toward  increasing  the  general  understanding  of  the  magneto- 
spheric  magnetic  field  and  increasing  the  gross  accuracy  of  the  models.  Efforts 
that  dwell  on  refining  some  of  the  details  of  present  mode's  might  profitably 
be  redirected. 

It  now  appears  to  be  well  established  that  magnetically,  the  magnetosphere  is 
open.  That  is,  a portion  of  the  earth's  magnetic  field  is  not  contained  in  the 
magnetosphere  but  passes  through  the  magnetopause  and  is  connected  to  inter- 
planetary magnetic  field.  It  therefore  would  not,  at  first  glance,  appear 
profitable  to  try  to  develop  closed  magnetic  field  models  (in  which  all  field 
lines  emanating  from  the  earth's  surface  return  to  the  earth's  surface  in  the 
opposite  hemisphere  without  going  outside  the  magnetopause.  However,  in  order 
to  accurately  determine  magnetic  field  properties  and  particle  behavior  near  the 
boundary  of  the  magnetosphere,  it  may  be  worthwhile  to  assume  a fixed  shape 
through  which  no  field  lines  can  penetrate  as  was  suggested  by  Roederer  (this 
technique  has  been  employed  by  Shabansky  and  his  colleagues).  An  alternative 
to  this  method  is  to  abandon  the  assumption  of  specular  reflection  off  the 
magnetopause  but  retain  the  pressure  balance  assumption  and  use  pressures 
obtained  from  a hydrodynamic  flow  model  such  as  developed  by  Sprei tor's  group 
from  wind  tunnel  tests  of  flow  past  an  objectshaped  to  represented  the  magneto- 
pause. Preliminary  work  on  this  problem  has  been  reported  by  Olson. 
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Electric  Fields 


The  current  situation  with  the  magnetospheric  electric  field,  £,  is  far  dif- 
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ferent  from  that  of  B.  Good  empirical  models  of  B and  the  source  currents  for 
-*■  -> 

B are  fairly  well  understood.  Progress  on  the  development  of  theoretical  B 

models  is  also  substantial.  However,  new  empirical  models  of  E are  just  now 

emerging  and  the  overall  picture  of  E is  still  cloudy  enough  that  no  good 

theoretical  models  seem  to  be  on  the  horizon.  There  are  several  reasons  for 

4 

this.  Throughout  the  magnetosphere,  E is  quite  variable  in  both  space  and 
time,  so  it  is  difficult  to  use  data  to  determine  average  characteristics  of 
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E.  Second,  most  work  on  E relies  heavily  on  its  indirect  measurement  using 
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particle  data  and  inferring  E from  the  V x B drift.  Third,  there  has  only 
been  limited  success  in  "mapping"  the  electric  field  from  one  region  of  the 
magnetosphere  to  another.  Typically,  measurement  of  the  scalar  potenfial 

-4 

associated  with  E is  determined  in  the  ionosphere.  The  assumption  tha..  field 
lines  are  equi potentials  is  then  made  in  order  to  map  this  electric  field  to  the 
equatorial  regions  of  the  magnetosphere.  There  are  at  least  three  possible 
sources  of  difficulty  associated  with  this  mapping  procedure.  First,  models  of 
the  magnetic  field  may  not  correctly  predict  the  observed  magnetic  field  line 

geometry.  Second,  if  induced  electric  fields  are  present,  the  assumption  fiat 
field  lines  are  equipotentials  is  no  longer  valid,  and  third,  anomalous  resis- 
tivity may  be  present  which  allows  electric  fields  to  have  a component  parallel 
to  the  direction  of  the  magnetic  field.  Cauffman  assessed  the  status  of  experi- 
ments in  operation  or  planned  for  the  direct  measurement  of  dc  electric  fields 
using  double  probes.  Many  long  boom-dyuble  probe  experiments  will  become  opera- 
tional during  the  IMS  era,  greatly  increasing  the  number  of  observations  of  t. 

However,  at  present,  most  dc  electric  field  observations  are  inferred  from  par- 
ticle data.  These  measurements  can  roughly  be  divided  into  those  made  in  the 
magnetosphere  and  in  the  ionosphere.  Higby  and  Hones  have  used  solar  particle 
motions  In  order  to  infer  magnetospheric  electric  fields,  while  Carpenter  and 
Rycroft  have  used  whistler  data  in  the  plasmasphere  to  infer  motions  that  are 
more  dynamo-like  than  convection- like.  In  the  ionosphere,  llellis  has  used 
Atmospheric  Explorer  data  to  develop  three-dimensional  particle  velocity 
measurements  and  from  them  high  latitude  convection  patterns.  The  Chatanika 
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data  was  discussed  in  some  detail  by  the  USCD  group.  Measurements  of  ion 
velocities  with  the  Cliatanika  radar  should  make  it  possible  to  determine  the 
role  of  discrete  auroral  arcs  in  the  convection  process  at  high  latitudes. 

In  a paper  read  by  Murayama  for  Mezeawa,  work  on  a reverse  current  flows  over 
the  polar  cap  was  reported.  This  suggests  motions  that  are  reverse  of  the 
standard  convection  patterns  inferred  at  high  latitudes.  Thus,  polar  cap  work 
and  plasmasphere  studies  have  suggested  deviations  from  the  classical  magneto- 
spheric  convection  patterns  and  indicate  once  more  the  high  degree  of  vari- 
ability in  the  magnetospheric  electric  field  and  also  possibly  some  physical 
processes  that  are  currently  not  well  understood. 

On  the  summary  panel.  Wolf  indicated  that  there  is  now  a fairly  complete  picture 
of  the  ionospheric  electric  field  and  it  may  be  possible  to  develop  an  empirical 
model  of  E in  the  ionospheric  regions.  This  model  might  be  similar  to  Mcllwain's 
electric  field  model  that  is  useful  in  equatorial  regions  of  the  magnetosphere 
and  would  have  as  inputs  to  it  such  variables  as  the  boundary  locations  of  the 
polar  cap.  Again,  the  usefulness  of  such  a model  would  depend  on  the  user  needs. 
Any  empirical  model  will  not  directly  help  to  further  our  understanding  of  mag- 
netospheric processes.  However,  deviations  between  such  a model  and  observations 
may  isolate  important  physical  processes  currently  unobserved.  Also,  such  a 
model  could  be  used  in  conjunction  with  available  magnetic  field  models  in  order 
to  study  charged  particle  phenomena  in  the  inner  magnetosphere  and  ionosphere. 

In  the  outer  magnetosphere,  the  electric  field  situation  is  much  more  complex 
both  from  a spatial  and  a temporal  point  of  view.  It  was  mentioned  before  that 
mapping  between  the  ionosphere  and  outer  magnetosphere  is  complicated  for  several 
reasons.  It  may  be  reasonable  to  try  to  couple  the  various  regions  of  the  magneto- 
sphere and  ionosphere  with  a semiempirical  model  in  which  magnetospheric,  iono- 
spheric and  plasmasphere  (whistler)  data  are  used  together.  There  was  a con- 
sensus, however,  that  there  is  no  chance  that  in  the  next  few  years  a theoretical 
model  of  E in  the  polar  cap  and  outer  magnetospheric  regions  will  be  developed. 

-*•  ■> 

Another  way  to  contrast  B and  E is  in  terms  of  their  sources.  Although  the 

source  currents  responsible  for  ft  are  well  known  (but  not  directly  measured)  to 
*♦ 

date  most  work  on  E has  looked  only  at  the  total  field  with  little  regard  for  its 
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sources.  Basically,  there  can  only  be  three  sources  for  E;  charge  separation, 
time  varying  magnetic  fields,  and  the  relativistic  transformation  of  the  solar 
wind  magnetic  and  electric  fields  into  a magnetospheric  reference  frame.  In 
order  to  develop  theoretical  models  of  the  magnetospheric  electric  field,  it 
will  be  necessary  to  understand  these  sources  of  E.  Of  these  three  source 
mechanisms,  probably  those  electric  fields  induced  by  time  varying  magnetic 
fields  are  most  amenable  to  understanding  and  quantitative  modelling.  Several 
sources  of  time  varying  magnetic  fields  in  the  magnetosphere  are  known;  changes 
in  the  standoff  distance  of  the  magnetopause  current  system,  substorm  related 
changes  in  the  tail  current  system,  storm  time  changes  in  the  ring  current,  and 
the  daily  variation  in  the  magnetic  field  produced  by  tilt  effects.  Hones  and 
Burgeson,  and  Birmingham  and  Northrup  did  early  work  in  general  problem  areas. 
Recent  procedures  for  calculating  induced  electric  fields  developed  by  Olson  and 
Pfitzer  have  shown  that  even  the  daily  wobble  of  the  earth's  magnetic  field  plus 
the  corresponding  wobble  in  the  magnetopause  currents  can  produce  a substantial 
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induced  electric  field.  The  same  sort  of  conments  made  on  B apply  to  E as 
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regards  user  needs.  The  models  of  E should  be  "global",  analytic,  etc. 

Emphasis  is  shifting  from  the  static  magnetosphere  to  the  study  of  its  dynamic 
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properties.  It  is  therefore  necessary  that  both  E and  B be  capable  of  reore- 

senting  observed  temporal  variations.  The  magnetospheric  substorm  is  an  example. 

In  order  to  properly  model  the  injection  of  charged  particles  during  substorms, 

it  is  necessary  to  employ  both  a time  varying  magnetic  field  model  and  a time 
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varying  electric  field  model.  Also,  during  time  varying  situations,  the  E and 

•4  -4 

B models  ere  coupled  in  that  the  time  variations  in  B will  induce  a component 
that  must  be  included  in  the  total  model  of  E.  Because  of  this  coupling  another 
requirement  must  be  made  on  all  future  magnetic  field  models  and  that  is,  in 
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addition  to  computing  B,  these  models  should  also  include  A,  the  vector  magnetic 
potential.  (The  induced  electric  field  is  best  computed  by  taking  the  time 

4 

derivative  of  A.)  The  need  for  having  a time  varying  electric  field  model  again 
points  out  the  importance  of  understanding  its  source  mechanisms  since  input 
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parameters  to  both  E and  B must  be  well  understood  if  the  models  are  to  generate 
meaningful  results. 


Charged  Particles 

In  the  session  on  charged  particles,  it  became  apparent  that  the  current  lack 
of  understanding  of  magnetospheric  boundary  phenomena  limits  the  use  of  field 
models  describing  charged  particle  behavior.  This  is  true  for  both  the  high 
and  low  energy  charged  particles.  The  access  of  low  energy  charged  particles 
is  not  well  understood,  nor  is  the  coupling  of  the  magnetospheric  magnetic  field 
with  the  magnetosheath  and  interplanetary  fields.  This  coupling  influences  the 
entry  of  energetic  particles  and  their  subsequent  motions  in  the  magnetosphere. 

It  is  possible,  however,  that  the  higher  energy  particles,  which  do  not  contri- 
bute  significantly  to  the  current  systems  that  produce  B can  be  used  to  study 
some  of  the  gross  properties  of  the  magnetospheric  magnetic  field. 

Masley  reported  work  on  comparisons  of  observed  solar  cosmic  ray  phenomena  and 
model  prediction.  Particle  trajectories  and  cutoffs  obtained  using  the  Olson- 
Pfitzer  model  were  in  good  agreement  in  the  noon  and  midnight  meridian  planes 
when  observations  were  used  for  the  same  range  of  Kp  values  as  the  model  is 
intended.  His  studies  of  cosmic  ray  entry  and  propagation  also  suggest  that 
energetic  particle  access  to  the  polar  cap  regions  is  strongly  influenced  by  the 
interplanetary  magnetic  field  direction.  Other  papers  on  high  energy  particles 
also  suggested  the  importance  of  accurate  magnetic  field  models  and  tests  that 
these  particles  can  offer  for  the  development  of  such  models.  It  was  suggested 
that  observations  of  high  energy  particle  fluxes  be  correlated  with  substorm 
time  to  help  develop  dynamic  models  of  the  magnetic  field.  Hoffman  and  Smith 
showed  that  particle  fluxes  observed  from  1 to  872  kilovolts  can  account  for 
most  of  the  depression  of  the  Inner  magnetosphere  observed  in  the  magnetic  field 
data  by  Sugiura  and  his  colleagues.  Buck  and  West  presented  detailed  pitch  angle 
distributions  from  which  they  derive  drift  shell  splitting  and  determine,  in- 
directly, the  shape  of  magnetic  field  lines.  Buck  has  examined  the  motions  of 
energetic  particles  near  noon  and  suggests  that  the  observed  pitch  angle  spectra 
can  be  explained  in  terms  of  a magnetic  field  that  has  minima  along  the  field 
lines  displaced  from  the  equatorial  region  as  much  as  60°.  Both  Konradi  and 
Mcllwain  suggested  that  the  region  between  6 and  10  earth  radii  where  particle 
injection  occurs  needs  more  detailed  study  in  order  that  it  be  properly  modelled. 
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It  is  in  this  region,  at  least  near  local  midnight,  where  the  magnetopause, 
quiet  time  ring  and  tail  current  systems  all  contribute  similar  amounts  to  the 
total  magnetic  field.  Murayama  discussed  energetic  electrons  that  were  observed 
near  the  high  latitude  dayside  magnetopause  several  hours  after  observed  sub- 
storm activity. 

All  of  these  particle  observations  must  eventually  be  accounted  for  by  improved 
field  models.  However,  it  appears  at  present  that  only  the  more  gross  features 
of  the  particle  and  field  phenomena  observed  in  the  magnetosphere  are  being  fit 
by  the  models.  These  detailed  magnetospheric  particle  features  will  provide 
tests  for  the  more  sophisticated  models  that  should  be  developed  in  the  near 
future. 

There  was  also  discussion  on  the  definition  of  the  last  closed  or  first  open 
field  line.  It  was  argued  that  this  was  best  done  using  charged  particle  data. 

A magnetometer  on  a polar  orbiting  satellite  cannot  by  Itself  distinguish  open 
from  closed  field  lines.  Magnetometer  data  can  best  be  used  to  define  average 
magnetic  field  models  and  then  use  the  model  geometry  to  define  the  boundary 
between  open  and  closed  lines.  However,  electrons  present  during  solar  events 
can  be  used  to  accurately  determine  the  position  of  the  first  open  field  line, 
although  this  technique  is  limited  because  of  the  Infrequent  occurrence  of 
solar  events,  It  nevertheless  can  be  used  to  determine  the  average  position 
(possibly  as  a function  of  Kp)  of  the  first  open  field  line. 

The  question  of  open  versus  closed  magnetosphere  as  regards  charged  particle 
entry  must  also  be  answered  in  the  affirmative.  That  is,  the  magnetosphere 
appears  to  be  always  "open"  to  the  entry  of  both  low  and  high  energy  particles 
over  some  portions  of  the  magnetopause.  This  is  true  even  if  the  magnetosphere 
Is  magnetically  closed.  This  Is  possible  because  some  particles  incident  upon 
the  magnetopause  enter  it  and  drift  through  the  magnetosphere  because  of  gradients 
in  the  magnetic  field  parallel  to  the  boundary.  For  lav  energy  particles  these' 
regions  Include  clefts  on  the  high  latitude  dayside  magnetopause  and  in  the 
equatorial  (plasma  sheet)  regions  of  the  tail.  Solar  and  galactic  cosmic 
rays  penetrate  all  regions  of  the  magnetopause.  It  was  concluded  that  the  magneto- 
sphere is  open,  both  in  the  magnetic  sense  and  from  the  standpoint  of  entry  of 
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charged  particles.  However,  present  models  do  not  explicitly  attempt  to  repre- 
sent this  open  nature  of  the  magnetosphere.  In  order  to  do  this  a much  better 
quantitative  understanding  of  the  magnetopause  must  be  obtained.  Such  a quanti- 
tative model  of  the  magnetopause  must  self-consistently  account  for  a (small?) 
component  of  magnetic  field  normal  to  the  boundary  and  a determination  of  which 
particles  are  deflected  from  tin’s  boundary  and  which  particles  are  able  to 
move  through  it. 

Conclusions 

The  goal  of  the  meeting  was  to  make  an  assessment  of  what  quantitative  models 
are  now  available  and  how  good  they  are,  to  determine  generally  what  data  sets 
are  available,  or  will  become  available  during  the  IMS  era,  and  then  to  deter- 
mine the  limitations  of  existing  models.  These  items  were  then  used  to  make 
recommendations  to  model  developers  and  model  users  for  work  that  should  be 
undertaken  during  the  IMS  period  to  extend  existing  models  and  to  develop  new 
ones.  A summary  of  the  present  situation  regarding  quantitative  models  is 
given  in  this  section  and  recommendations  for  work  in  the  near  future  given  in 
the  next. 

Modelling  efforts  in  magnetic  fields,  electric  fields,  and  charged  particles 
were  described.  Clearly,  efforts  to  model  the  magnetospheric  magneti c field 
has  made  the  most  progress.  However,  it  is  clear  that  an  "international" 
reference  model"  is  not  yet  at  hand.  There  are  several  models  available  that 
do  specific  things  well,  but  no  single  model  that  meets  all  of  the  user  require- 
ments. The  current  models  range  from  physical  to  entirely  empirical.  The 
utility  of  a given  model  depends  on  the  needs  of  the  user.  This  was  one  of  the 
key  problem  areas  isolated  at  the  meeting.  Although  different  users  typically 
would  ask  for  the  same  output  from  a model,  the  inputs  might  vary  considerably. 
For  example,  a user  who  is  trying  to  understand  a particular  magnetospheric 
phenomenon  may  only  require  that  a model  input  parameter  becomes  available 
several  weeks,  or  even  months,  after  the  time  in  question  (for  example,  Kp 
or  many  of  the  other  indicies  of  magnetic  activity)  while  there  are  many 
"more  practical"  users  who  require  input  parameters  which  can  themselves  be 
monitored  in  real  time  (for  example,  the  solar  wind  dynamic  pressure,  the 
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direction  of  the  interplanetary  magnetic  field).  Thus,  determining  the  proper 
input  parameters  for  a model  is  not  a trivial  problem  is  the  needs  of  these 
several  users  are  to  be  met.  It  is  necessary  that  general  input  parameters 
for  the  models  be  developed.  This  means  that  in  addition  to  the  location  in 
question,  inputs  to  a model  that  will  give  back  the  vector  magnetic  field  m#st 
also  include  parameters  like  the  standoff  distance  of  the  magnetopause  which 
can  be  related  either  to  l(p  (after  the  fact)  or  in  real  time  to  solar  wind 
parameters.  Both  model  and  data  user  communities  indicated  that  a data  sam- 
pling rate  of  1 per  minute  suffices  for  most  tests  and  applications.  Some 
users,  of  course,  use  much  lower  sampling  rates.  Po wer  suggested  a rate  of 
1 per  hour  for  some  practical  applications. 


Many  indices  and  real  time  input  parameters  were  discussed.  They  included  the 
location  of  the  auroral  oval  and  the  electron  trough,  the  size  of  the  polar  cap 
the  Q index  Kp,  A^-,  D^,  the  vector  specification  of  the  interplanetary  magneti 
field,  the  various  solar  wind  parameters,  and  the  tilt  angle.  For  input  para- 
meters to  the  models  to  have  general  utility,  it  is  necessary  that  they  relate 
both  to  parameters  that  are  readily  available  (such  as  Kp)  but  not  necessarily 
in  real  time,  and  to  other  environmental  features  that  can  be  monitored  in  real 
time  (such  as  the  dynamic  solar  wind  pressure).  Thus  the  standoff  distance 
(which  gives  an  indication  of  the  size  of  the  magnetosphere  and  concurrently 
the  strength  of  the  magnetopause  current  systems)  might  be  an  appropriate  model 
input  parameter.  It  would  be  then  left  to  the  user  to  relate  this  particular 
model  input  parameter  to  another  quantity  convenient  for  his  particular  use 
(for  example,  Kp,  or  the  solar  wind  dynamic  pressure).  Many  such  input 


parameters  must  be  decided  upon.  This  must  be  done  v/ith  both  user  needs  and 
data  availability  kept  in  mind. 

Although  there  were  several  complaints  about  ava. .ability  of  good  data  for  use 
in  the  construction  of  empirical  models,  and  as  a monitored  input  to  some 
models,  it  is  apparent  that  the  modelling  community  must  make  the  best  of  what 
is  available.  For  example,  it  was  suggested  that  vector  aB  measurements  (aB 
is  the  difference  between  the  observed  magnetic  field  and  that  portion  produced 
by  the  earth's  magnetic  field  as  given  by  an  accurate  model  of  the  main  field) 
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be  used  instead  of  scalar  aB  measurements.  When  vector  aB  measurements  become 
generally  available,  they  will  be  used  as  inputs  to  models  and  also  to  test 
models.  To  date,  however,  the  scalar  aB  measurements  have  been  used  effectively 
to  understand  current  systems  in  the  inner  magnetosphere  and  to  develop  quanti- 
tative models  of  the  inner  magnetospheric  magnetic  field. 

4- 

The  general  availability  of  data  and  the  fact  that  B can  be  measured  directly 
has  served  to  expedite  efforts  to  quantitatively  model  magnetospheric  magnetic 
fields.  Efforts  to  model  the  electric  field,  on  the  other  hand,  have  suffered 

4 

from  a lack  of  direct  measurement  of  E and  its  extreme  variability,  both  in  time 
and  space.  Efforts  to  "map"  the  electric  field  between  different  regions  of  the 
magnetosphere  and  thus  compare  simultaneous  data  sets,  have  met  with  only  partial 
success,  due,  in  part,  to  this  variability.  Hopefully,  as  several  double  probe 
experiments  become  available  to  directly  measure  the  dc  electric  field,  some  of 
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these  problems  will  be  resolved,  and  the  effort  to  model  E in  terms  of  its 
various  sources  can  be  accelerated. 

For  the  future  development  of  models,  the  various  users  should  be  kept  in  mind. 
Users  vary  from  those  interested  only  in  an  understanding  of  the  magnetosphere 
(and  an  interest  only  in  completely  physical  theoretical  models)  to  those  prac- 
tical users  who  are  interested  in  models  simply  as  "black  boxes"  that  yield 
quantitatively  accurate  information.  Between  these  two  extremes  there  is  an 
almost  continuous  spectrum  of  users.  These  include,  for  example,  the  experimenter 
who  needs  to  know  the  location  of  the  foot  of  the  magnetic  field  line  that  passes 
through  his  satellite,  and  some  monitoring  activities  that  can  be  done  within 
several  hours  of  real  time  observations.  Some  of  the  practical  uses  of 
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magnetospheric  models  include  the  understanding  and  prediction  of  spacecraft 
cfiarging  and  arcing,  coupling  of  the  magnetosphere  to  ionosphere  arid  associated 
radio  communications  problems,  the  coupling  of  the  magnetosphere  and  ionosphere 
to  the  upper  atmosphere  and  the  influence  of  the  upper  atmosphere  on  satellite 
trajectories  and  lifetimes,  and  the  effects  of  trapped  radiation  and  radiation 
present  during  solar  cosmic  ray  events  on  various  space  hardware.  Hildebrand 
suggested  that  the  scientific  community  keep  in  mind  the  needs  of  the  practical 
users  because  they  can  exert  considerable  influence  on  the  determination  of 
funding  levels  in  several  government  agencies. 

Recommendations 

In  order  to  best  meet  current  and  future  user  needs  for  quantitative  models  of 
magnetospheric  phenomena,  several  recomnendations  are  made  here  that  should  have 
special  significance  during  the  IMS  era.  During  the  IMS,  increased  attention 
will  be  focused  on  all  aspects  of  magnetospheric  phenomena.  These  efforts  will 
include  sets  of  satellite  and  ground  based  experiments  and  monitoring  stations 
the  procurement  of  many  observational  data  sets  and  the  development  of  both 
physical  and  empirical  quantitative  models.  The  first,  and  probably  most 
general  recommendation,  is  that  a model  situation  center  be  formed  - in  the 
same  manner  as  the  experiment  situation  center.  Such  a center  would  serve  two 
important  purposes.  First,  it  would  permit  a means  for  closely  monitoring 
progress  in  the  development  of  models  and  a means  for  evolving  new  sets  of 
modelling  recommendations  as  better  data  sets  become  available  and  more 
sophisticated  quantitative  models  are  developed.  Second,  such  a center  would 
provide  a means  for  interaction  between  the  magnetospheric  physics  community 
and  the  many  users  who  need  to  know  quantitatively  the  effects  of  the  magneto- 
sphere (possibly  coupled  with  the  ionosphere  and/or  atmosphere)  on  their  hard- 
ware systems.  It  should  be  remembered  that  these  users  can  provide  much  needed 
incentive  to  various  agencies  to  continue  their  support  of  magnetospheric 
research.  In  the  eyes  of  these  users,  the  quantitative  model  is  the  end  result 
of  the  magnetospheric  research.  The  model  is  employed  by  the  user  to  cost 
effectively  find  solutions  to  environmental  effects  on  his  hardware  systems. 

A modelling  situation  center  would  also  provide  an  interface  between  modelling 
and  observational  data  collection  efforts. 


* 

* 


More  specific  recommendations  concerning  the  future  development  of  quantitative 
magnetospheric  magnetic  field  models  were  also  made.  Several  of  these  recom- 
mendations  apply  as  well  to  the  development  of  electric  field  models  and  some 
even  to  the  development  of  charged  particle  models.  These  recommendations  can 
be  roughly  divided  into  two  categories;  requirements  that  future  models  should 
meet  and  quantitative  tests  that  they  should  satisfy.  General  requirements 
that  any  future  quantitative  magnetospheric  magnetic  field  model  should  satisfy 
include: 

1.  It  should  represent  the  magnetic  field  contributions  from  all  important 
magnetospheric  current  systems  including  those  currents  distributed 
throughout  the  inner  magnetosphere  (the  quiet  and  disturbed  ring 
currents)  and  in  the  plasma  sheet  region  of  the  tail  (the  tail  current 
system) . 

2.  The  dependence  of  each  of  these  current  systems  on  the  tilt  angle 
(the  complement  of  the  angle  between  the  solar  wind  direction  and  the 
geomagnetic  dipole  axis)  should  be  included. 

3.  Variability  In  each  of  the  current  systems  should  be  taken  into  account 
(e.g.,  the  shape  and  strength  of  the  magnetopause  current  systems  as  a 
function  of  standoff  distance). 

4.  The  representation  of  B should  be  analytic. 

-+  -+ 

5.  The  magnetic  vector  potential  A should  be  computed  together  with  B. 

6.  The  model  should  be  valid  over  a large  portion  of  the  observed  magnetosphere. 
In  addition,  there  are  some  specific  user  requirements: 

7.  The  amount  of  computer  time  needed  to  represent  3 should  be  minimized. 

8.  Input  parameters  to  the  computer  subroutine  should  be  easily  obtainable 
from  user  chosen  indices  or  monitored  parameters. 
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The  first  three  requirements  would  satisfy  most  of  the  complaints  made  against 

existing  models:  they  do  not  include  the  effects  of  all  current  systems  and/or 

tilt  effects,  and  most  important,  they  do  not  do  a good  job  of  describing 

dynamic  magnetospheri c behavior.  A dynamic  model  must  of  necessity  have  input 

parameters  that  are  related  to  tne  variability  of  the  several  current  systems 
-► 

contributing  to  11.  In  order  to  use  the  model  to  describe  charged  particle 
behavior,  it  is  necessary  that  it  be  both  analytic  and  differentiable.  Such 
a model  of  B should  then  yield  both  the  motions  of  charged  particles  (ifforn 

-*•  -►  f 

vB)  and  currents  distributed  throughout  the  magnetosphere  (from  vxB).  It  is 

necessary  to  compute  the  vector  potential  A in  order  to  determine  the  electric 

-*■  -*  *+ 
field  induced  by  temporal  variations  in  B.  A can  also  be  used  to  determine  B 

■+ 

if  it  is  demanded  that  v*B  be  precisely  0.  A computer  representation  of  li 
must  be  extremely  fast  if  it  is  to  be  used  to  compute  magnetic  field  lines, 
low  energy  particle  motions  and  drift  shells,  and  the  precise  trajectories 
of  high  energy  particles. 

In  order  to  assess  the  quantitative  accuracy  of  any  model  meeting  the  above 
criteria,  it  is  necessary  to  have  available  quantitative  tests.  It  is  recom- 
mended  that  a standard  set  of  tests  be  applied  to  all  B models  (it  is  under- 
stood that  some  specialized  models  would  not  do  as  well  as  others  in  satisfying 
these  tests  but  still  be  the  most  useful  for  certain  specific  user  needs).  Several 
existing  models  have  been  tested  by  comparing  their  output  with  quantitative 
observational  data  of  known  magnetospheri c phenomena.  These  data  include: 

1.  Magnetic  field  observations  in  different  regions  of  the  magnetosphere 
(scalar  <\B  plots,  the  decay  of  the  tail  field  in  its  equatorial  and  lobe 
regions,  the  field  geometry  at  high  latitudes  - especially  near  tlje 
magnetopause,  the  latitude  of  the  last  closed  magnetic  field  line). 

2.  The  daily  variations  in  the  synchronous  magnetic  field  during  quiet 
and  disturbed  periods. 

3.  The  trapping  boundary  of  adiabatic  particles. 

4.  The  low  latitude  cutoff  boundary  for  solar  cosmic  rays. 
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5.  Field  line  shape  as  determined  from  barium  releases. 

Other  observational  data  sets  that  should  be  used  to  test  models  in  the  near 
future  include: 

6.  Vector  aB  measurements 

7.  Detailed  pitch  angle  information  on  trapped  particle  phenomena. 

Several  of  the  same  requirements  can  be  applied  to  the  electric  field  modelling 
efforts.  The  greatest  problem  with  the  electric  field  to  date  appears  to  be 
the  lack  of  data  in  the  inner  and  outer  regions  of  the  magnetosphere.  However, 
because  there  is  currently  a fair  amount  of  data  on  the  ionospheric  electric 
field,  it  is  recommended  that  an  empirical  model  of  the  ionospheric  electric 
field  be  developed.  Such  a model  would  at  least  indirectly  help  the  under- 
standing of  electric  fields  throughout  the  magnetosphere.  The  problem  of  data 
in  the  inner  and  outer  magnetosphere  regions  should  largely  be  overcome  by  the 
deployment  of  several  long  boom,  double  probe  experiments  to  measure  the  dc 
electric  field.  As  this  data  becomes  available,  it  is  recommended  that  the 
various  sources  of  l be  kept  in  mind  in  trying  to  explain  and  model  the  magneto- 
spheric  electric  field.  For  example,  magnetospheric  electric  field  observations 
could  be  checked  for  the  existence  of  a longitudinal  component  which  is  dependent 
on  universal  time.  Such  a component  in  the  total  field  could  be  induced  b>  the 
daily  wobble  of  the  earth's  dipole  field.  The  availability  of  simultaneous  data 
from  different  regions  of  the  magnetosphere  and  ionosphere  will  also  shed  light 
on  the  mapping  problem. 

» 

It  is  probably  the  lack  of  understanding  of  the  various  boundaries  in  the  magneto 
sphere  and  ionosphere  that  causes  the  largest  problem  for  the  physical  modelling 
of  electric  fields.  Thus  the  proper  determination  of  the  properties  of  tike 
magnetopause  again  appears  important  (It  was  mentioned  above  that  the  magneto- 
pause properties  must  be  understood  in  detail  in  order  to  understand  the  coupMng 
of  the  interplanetary  magnetic  field  with  the  geomagnetic  field  and  in  order  to 
properlv  understand  and  model  the  entry  of  charged  particles  into  the  magneto- 
sphere). Thus  i t is  recomnended  that  some  attention  be_  directed  toward  the 
!nJ*SlP.6 top a us  e , such  that  ^lectHc_  and  magnetic  fields  and  particles  jQiienoincna 
both  on  .and  through  that  boundary  cap  be  quant A tati vely ’ unders  tqpd. 
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There  are  now  many  data  sets  that  detail  charged  particle  behavior  in  the 
magnetosphere.  It  would  be  of  general  help  to  the  electric  and  magnetic 
field  contingent  of  the  magnetospheri c community  if  models  were  to  be 
developed  of. several  of  the  gross  particle  features  in  the  magne*  sphere 
(e.g.,  the  properties  of  the  plasma  sheet  and  the  plasmasphere).  It  is 
probably  necessary  to  have  a first  order  set  of  models  (of  the  electric,  and 
magnetic  field,  and  of  charged  particle  structures)  in  order  to  finally  consider 
their  mutual  interaction  and  begin  the  development  of  a true  quantitative  mag- 
netospheric  model.  It  might  be  argued  that  models  of  the  magnetic  field,  or 
of  charged  particle  phenomena,  are  merely  different  ways  of  representing  the 
same  magnetospheri c phenomena.  A quantitative  magnetospheri c model,  then, 
may  at  least  be  defined  as  one  that  self-consistently  and  simultaneously  con- 
siders the  magnetospheri c magnetic  and  electric  fields  and  charged  particle 
phenomena. 

It  was  also  recommended  that  considerable  effort  be  devoted  toward  the  deter- 
mination of  proper  in p_ut_ parameters  to  all  of  the  models . This  is  i mp e r a t i v e 
if  the  models  are  to  meet  the  various  user  needs.  These  input  parameters  must 
be  chosen  so  as  to  make  the  models  dynamic.  In  the  case  of  U,  input  parameters 
are  needed  to  describe  each  of  the  various  current  systems.  It  may  be  that  for 
some  user  applications  existing  indices  are  adequate  (e.g.,  as  an  input  to 
described  ring  current  strengths).  However,  it  appears  that  no  parameters 
presently  are  adequate  to  describe  the  dynamics  of  the  tail  current  during 
magnetospheri c subs  forms.  It  was  suggested  that  the  size  of  the  polar  can 
might  provide  a good  indication  of  the  strength  and  location  of  the  tail 
currents.  It  is  clear  that  a given  model  "can  be  no  better  than  its  input 
parameters".  Thus  the  selection  and  use  of  the  proper  indices  and/or  monitored 
data  for  input  parameters  must  be  considered  an  important  part  of  the  develop- 
ment of  quantitative  models  of  magnetospheri c phenomena. 
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The  thermospheric  model  In  our  August  1974  report  to  AFOSR  [Olson  et  al.,  1974] 
Includes  a quantity,  f(t),  which  represents  the  diurnal  density  variation  at 
the  equator.  In  Table  1 of  that  report,  f(t)  Is  represented  by  low-latitude 
data  collected  by  0G0-6  [Hedln,  at  al.,  1974].  This  table  has  now  been  compared 
with  measurements  obtained  by  accelerometers  on  the  two  equatorial  satellites, 

San  Marco  II  [Brogllo,  ’971]  and  III,  [Brogllo  et  al.,  1974],  and  with  the 
Paetzold  [1963]  model,  which  was  based  on  data  from  seven  early  satellites.  All 
four  data  sets  were  compared  at  an  altitude  of  450  km  (In  the  case  of  the  two 
« San  Marcos,  this  required  a moderate  amount  of  extrapolation).  The  four  data 

sets  have  been  averaged  to  obtain  the  values  plotted  as  open  circles  In  Figure  1. 

By  using  Karl  Pfltzer’s  universal  fitting  program,  ORTHON,  the  following  Fourier 
series  has  been  fitted  to  the  points: 

* 

f(t)  « 0.994  + 0.545  cos  C(t  - 14.745)  ♦ 0.102  cos  2C(t  - 1.838)  f 

+ 0.0154  cos  3C(t  - 23.098)  - 0.00667  cos  4Ct+ 0.0161  cos  5C(t  - 0.777) 
where  C • ir/12  and  t Is  local  time  In  hours. 

The  equation  Is  also  shown  In  Figure  1. 
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In  addition  to  the  seasonal  variation,  there  is  a semiannual  variation  in 
atmospheric  density  which  was  discovered  by  Paetzold  and  Zschorner  [l9oij. 
Its  existence  was  at  first  denied  by  Jacchia  and  some  other  workers,  but 
it  is  now  firmly  established.  The  semiannual  variation  has  been  detected 
from  90  km  well  into  the  exosphere.  Some  measurements  of  the  ratio  of  the 
October  maximum  to  the  July  minimum  are  shown  in  Figure  1. 


Z (KM) 


Figure  1 . Ratio  of  Density  at  October  Maximum  to  July  Minimum.  Plotted  as  a Function  of  Altitude 
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Jacchia's  [1971]  analytic  approximation,  which  is  based  on  data  from  sever, 
satellites  and  Cook's  rocket  data,  is  shown  dashed.  Measurements  between  90 
and  215  km,  which  were  largely  lacking  from  Jac<_hia's  data  sample  are  shown; 
also  the  range  of  data  near  1*00  and  1100  km  for  a number  of  years  reported 
by  Cook  [1970],  and  three-year  averages  at  900  and  1070  km  published  by  Cook 
[1972].  Jacchia's  formula  does  not  fit  th'  daJ a below  200  km.  Sven  at 
higher  altitudes,  the  variability  of  the  tnef.:cx  -d  values  does  not  justify 
the  elaborate  formula  constructed  by  him.  A ampler  formula,  represented 
the  solid  curve,  has  therefore  been  constructed  using  K.  A.  Pfitzer's  rout 
ORTHOH. 

In  the  revised  formula,  R(Z)  represents  the  ratio  of  the  October  density  maxi- 
mum to  the  July  minimum!. 

R(Z)  = 0.98  + 0.27  x 10*2  Z - 0.85  x 10‘6  Z2 
-0.59  x 10‘9  Z3. 


The  tine  dependence  is  approximated  by  a sum  of  cosines  of  period  1 year, 
1/2  year,  1/3  year,  etc.  The  phases  of  the  cosine  terms  are  adopted  from 
/olland  et  al.  [1972].  The  result  plotted  in  Figure  2,  using  ORTHOH,  is: 

G(D)  = O.1U3  cos  (K( D— U ) ) + 0.239  cos  (2K(D-109)) 


+ 0.0U1*  cos  (3K(D-66)) 


where  K ^ 


and  D is  the  day  number  in  the  year. 


Aduing  a fourth  harmonic  term  does  not  significantly  improve  the  fit.  The 
semiannual  variation  of  density  is  then  represented  by  multiplying  the 
expression  for  density  by  the  factor 

[1  + (R(Z)  - 1)  G(D) ]. 
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Figure  2.  Time  Dependence  of  the  Density  When  the  Amplitude  Ratio,  R * 2. 
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